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Summary 

Reaction rate coefficients and thermodynamic and transport properties are pro- 
vided for the 11-species air model which can be used for analyzing flows in chemical 
and thermal nonequilibrium. Such flows will likely occur around currently planned 
and future hypersonic vehicles. Guidelines for determining the state of the surround- 
ing environment are provided. Approximate and more exact formulas are provided 
for computing the properties of partially ionized air mixtures in such environments. 

Introduction 

Currently envisaged transatmospheric and aeroassist missions (refs. 1,2,3, and 4) 
have created a resurgence of interest in the aerothermodynamic design of hypersonic 
vehicles. However, the velocities and altitudes at which these proposed craft would 
operate are different, and sometimes more severe, than have been experienced in the 
past. As a result, the nonequilibrium flow environment which will surround these 
vehicles will considerably impact the vehicle aerodynamics, thermal loads, and propul- 
sion system efficiency. Since such an environment is difficult to simulate in current 
ground based-test facilities, the design of these future vehicles will rely heavily on 
numeric calculations. In turn, these calculations will require a good understanding of 
the physical modelling required to simulate these phenomena. 

Under hypersonic flight conditions, a vehicle travelling through the atmosphere 
will excite the air which flows around the body to very high temperatures as kinetic 
energy from the vehicle is transferred. Depending on the flight velocity, various 
chemical reactions will be produced behind a shock wave as shown in figure 1 (which 
is adapted from ref. 5) for the stagnation region of a 30.5 cm (1 ft) radius sphere. 
These reactions will affect the properties of air and cause considerable deviation from 
those of a thermally and calorilicaily perfect gas. A vehicle flying through the higher 
reaches of the atmosphere at high velocities may also experience thermal nonequili- 
brium (fig. 1), since the lower density reduces the collision frequency, and the high 
velocity results in smaller transit times for the air molecules. Both of these processes 
create a delay in the equilibration of translational, rotational and vibrational modes of 
the thermal energy. Under these conditions, the modelling of the air chemistry 
requires a multi-temperature approach in contrast to classical single temperature for- 
mulations. 

Figure 1 delineates four regions (I through IV), showing when the various chem- 
ical activities are initiated at a given altitude and velocity. Similarly, it also depicts 
through Regions A, B and C the initiation of chemical and thermal nonequilibrium 
processes for different velocity and altitude conditions. This figure clearly shows that 
the reaction rate coefficients and thermodynamic and transport properties would 
change continuously for a given flight trajectory. For example, in regions A and B 
(i.e., before initiation of thermal nonequilibrium), the specific heat at constant pres- 
sure C p would change as follows: 

C p = constant in Region I, before the excitation of vibrational energy mode 

= C p (T ) in Region I, after the excitation of vibrational energy mode and 
before the dissociation of oxygen 
= Cp(C t ,T) after the dissociation of oxygen 

Similarly, the equation of state would change along the flight trajectory as the 
thermal equilibrium and thermal nonequilibrium regions (ref. 6) are traversed. 


I 
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p = p — zr—r in Regions A and B 
M 

R 

= p — ^-T lr in Region C 

M 

In numerical simulations, the thermodynamic and transport properties and reac- 
tion rate coefficients (in the case of finite rate chemistry) arc typically required. It is 
obvious from the previous discussion that these properties and the equation of state 
should be evaluated carefully when chemical and thermal nonequilibrium conditions 
exist in the flowfield around a hypersonic vehicle. Under chemical and thermal equili- 
brium conditions, the transport and thermodynamic properties of high temperature air 
and its components are well documented in the literature (refs. 7,8, and 9). However, 
for flows with finite-rate chemistry the individual species properties and appropriate 
mixing laws which are required are not as well established. For example, in a partially 
ionized gas mixture, the conventional mixing laws (refs. 10 and ll) developed for 
non-ionized mixtures can not be extended to higher temperatures without consider- 
able error (ref. 12). 

The purpose of this report is to provide thermodynamic and transport properties 
and the reaction rate coefficients of the most important reactions for the 11 consti- 
tuent species of air (N ,0 ,N 2 ,0 2 ,N0 ,A r + ,0 + ,A r 2 . 0 } ,NO + ,e~) for temperatures up to 
30000 K . Approximate and more exact mixing laws are also provided for partially 
ionized gas mixtures. Sources of the input data used in the calculation of various 
flowfield properties are identified. Appropriate formulas are provided for using these 
properties in computations of flows with thermal nonequilibrium. 


Symbols 
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b ,r 


A 




coefficients of polynomial curve-fits for thermodynamic proper- 
ties, n= 1,2. ...7 

coefficient in the Arrhenius form of backward reaction rate con- 
stant 

coefficient in the Arrhenius form of forward reaction rate con- 
stant 


A U 


A 'j 

a d„’ b d„' C D,j' d 5 „ 


nondiagonal matrix elements of the first Chapman-Enskog for- 
mula 

ratio of collision cross sections, ■j [A) 

curve-fit coefficients for diffusion coefficient D tJ 


a k, b k, C Kf p Kf ,E Kf . 

coefficients of polynomial curve fits for frozen thermal conduc- 
tivity of species i 

A ^ . , B ^ , , C . , E ^ . 

curve-fit coefficients for viscosity of species i 


A pro,!), B fru.i), Crru.i), Dq(1,d 

w <> >i ■> n 

curve-fit coefficients for collision cross section ' 
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A jyc.2.2), B ^-(2.2), C 2.2), D -ft (1.1) 


A a •, 5 « •, C R • , Z) 

B 'iJ 

B b ,r 

B f.r 

C, 


'Pi 


'Pf 


'p,t 


D u 

D U 

Eb.r 

E f.r 


F° 

h 

hi 

h >, 


vib 


(A h/) T 


rtf 


J- 

K 

K d 

K e 

«: 


K 


el 


B’j 


curve-fit coefficients for collision cross section 0^ 2,2) 
curve-fit coefficients for collision cross section ratio fi ; * 

ratio of collision cross sections, (5iT^ l,2) - 4£7 1 $ 1,3 Vi2y 1,1) 
temperature exponent for backward reaction rate constant 
temperature exponent for forward reaction rate constant 

mass fraction of species i 

r 

dh 


specific heat at constant pressure, 


ar 


, cal! gm -mole- K 


diffusion specific heat at constant pressure, 

ns ac ( - 




, call gm- mole- K 


p 


frozen specific heat, 
NS 


Cp f ~ X C iCp,i 


i* l 


ii. v A i£i 
ar -tr ar 


, call gm- mole- K 


dhj 

IF 


iP 


, call gm-mole-K 


specific heat of species i, 

coefficient of self-diffusion, cm 2 / sec 
binary diffusion coefficient, cm 2 / sec 

activation energy for the backward reaction r, erg I gm- mole 
activation energy for the forward reaction r, erg/ gm-mole 

free energy of species i at 1 atm pressure (standard state), 

call gm-mole 

NS 

enthalpy of mixture, cal/ gm-mole 

i*i 

enthalpy of species i, cal/ gm-mole 

vibrational component of enthalpy of a molecular species i, 
cal/ gm-mole 

standard heat of formation of species i at temperature T ref , 
cal! gm-mole 

diffusion mass flux of species i, gm I cm 2 - sec 

total effective thermal conductivity of mixture in therm o- 
dynam ic equilibrium, Kf + K d or Kf + K r , cal/cm-sec-K 

diffusion component of thermal conductivity of mixture defined 
by equation (A 14), cal/cm-sec-K 

thermal conductivity of electrons, cal/cm-sec-K 

thermal conductivity of electrons due to electron-electron colli- 
sions only, cal/cm-sec-K 

component of thermal conductivity of mixture due to electron 
excitation, cal/cm-sec-K 
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K f 

K fJ 


Kim.i 

K r 

K,r*l X) 

K 


K 


tr,i 


K 


vib 


k 

^ b ,r 

k f.r 

Le f.i 

Le ri 

Mi 

M 

n 

Pr r 

P 


Pvt b 


/?, 



frozen thermal conductivity of mixture in thermodynamic equili- 
brium, Ktr + K^t, call cm- sec- K 

frozen thermal conductivity of species i in thermodynamic 
equilibrium, K, ri + K int i , call cm-sec-K 

internal component of frozen thermal conductivity of a mixture 
in thermodynamic equilibrium, K rot + K vib + K el , call cm-sec-K 

internal component of the frozen thermal conductivity of species 
i in thermodynamic equilibrium, call cm-sec-K 

reaction component of thermal conductivity of mixture defined 
by equation (A16), call cm-sec-K 

translational thermal conductivity of mixture from first 
Chapman-Enskog approximation, call cm-sec-K 

translational thermal conductivity of mixture without contribu- 
tions due to electron-heavy particle collisions, K lr - K e , 

call cm-sec-K 

translational component of thermal conductivity of species i, 
call cm- sec- K 

vibrational component of thermal conductivity of mixture, 
call cm- sec- K 

Boltzmann’s constant, 1.38xl0 -16 erg/K 

backward reaction rate coefficient for reaction r, cm i l mole-sec 
or cm 6 7 mole 2 - sec 

forward reaction rate coefficient for reaction r, cm* I mole- sec 

frozen binary Lewis number, pC Pf DijlKf 

reactive binary Lewis number, pC p D t jlK 

molecular weight of species i, gm/gm-mole 

molecular weight of mixture, gm/gm-mole 

Avogadro’s number, 6.0225x 10 23 molecules / gm - mole 

number density, particles / cm* 

frozen Prandtl number, C P/ \i/Kf 

reactive Prandtl number, C p \il K 

pressure, dynelcm 2 

£th component of the overall heat flux vector, call cm 2 -sec 

k th component of the electron heat flux vector due to electron- 
electron collisions, call cm 2 -sec 

&th component of the vibrational heat flux vector due to 
molecule-molecule collisions, call cm 2 - sec 

universal gas constant, 1.987 call gm-mole-K 
temperature under thermodynamic equilibrium, K 
characteristic reaction temperature for the backward reaction r, 
E b<r lk,K 
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T e 

T et 


T 

T 

T 

T 

T 


ref 

rot 

tr 

ve 

vib 


V f 

VVj 

Xi 


Xi 


O-i.r 

hr 

yj 

n (,) .n 



rrp> 


fi-f.w ir.< u) 

ij ij 


characteristic reaction temperature for the forward reaction r, 
Ef tr lk, K 

electron temperature, K 
electronic excitation temperature, K 
reference temperature, 298.15 K 
rotational temperature, K 
translational temperature, K 

vibrational-electron-electronic excitation temperature, K 
vibrational temperature, K 
time, sec 

£th component of diffusion velocity of species i, cm I sec 
mass rate of formation of species i, g/ cm 3 - sec 
concentration of species i, moles I volume 
k th component of general orthogonal coordinate system 
mole-fraction of species i 
stoichiometric coefficients for reactants 
stoichiometric coefficients for products 
mole mass ratio of species j, Xjf p 

viscosity of mixture from first Chapman-Enskog approximation, 
gm! cm -sec 

viscosity of species i, gm I cm- sec 
density of the mixture, gmtem 3 
density of species i, gm/cm 3 

average collision cross section (for diffusion, viscosity, and 
translational, internal and reaction components of thermal con- 
ductivity) for collisions between the species i and /, A 2 ; 

1 A= 10- 8 cm 


average collision cross section (for viscosity and translational 
component of thermal conductivity) for collisions between the 
species i and j, A 2 


average collision cross sections (for translational component of 
thermal conductivity) for collisions between the species i and /, 

A 2 


Subscripts: 

b 


e 

el 

f 

i 


backward reaction 
electron 

electronic excitation 
forward reaction 
species / 
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J 

P 

r 

ref 

rot 

tr 

ve 

v 


species j 

constant pressure 

rth reaction 

reference condition 

rotational energy mode 

translational energy mode 

vibrational-electron-electronic energy mode 

total number of species, a NS 


Abbreviations: 

AOTV 

NASP 

NIR 

NJ 

NR 

NS 


Aeroassisted Orbital Transfer Vehicle 

National Aero-Space Plane 

total number of independent reaction 

sum of the reacting species (NS) plus the number of catalytic 
bodies 

total number of reactions 
total number of species 


Chemical Kinetic Model and Reaction Rates 


When chemical reactions proceed at a finite rate, the rate of production terms 
appear in the energy equation when formulated in terms of temperature and in the 
species continuity equations (refs. 13,14,15). For a multicomponent gas with NS 
reacting chemical species and NR chemical reactions, the stoichiometric relations for 
the overall change from reactants to products are: 

NJ NJ 

IP i,r*i (1) 

1 y <-l 

K b,r 


where r=l,2,...,NR and NJ is equal to the sum of the reacting species (NS) plus the 
number of catalytic bodies. The quantities a,,. and (3, r are the stoichiometric 
coefficients for reactants and products, respectively, and k f r and k b r are the forward 
and backward rate constants. The quantities X t denote the concentrations of the 
chemical species and catalytic bodies in moles per unit volume. The catalytic bodies 
(NJ-NS) may be chemical species or linear combinations of species that do not 
undergo a chemical change during the reaction. 


The net mass rate of production of the ith species per unit volume resulting from 
all the reactions NR may be obtained (ref. 14) from 


w« 


NR 


r- 1 


r 


dX, 




dpi_ 
dt r 

l Jr 


(2a) 


or 


w, 


NR 


M,X(P i,r-a.i.rKRf.-R b ,r) 

r= 1 


(2b) 
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where 


NJ 


Rf, = k/.rU(y ,P) " 

7- i 

NJ „ 

= **.rII(Y;P) 

7-1 


Here, the mole-mass ratio y ; is defined as 

X 


j= 1.2.....NS 


_j_ =£l 
P Mj 

NS 

(j-NS),i Y» j-NS+l,...,NJ 

i=i 


(2c) 

(2d) 


(2e) 


The constants Z (j-NS).i are determined from linear dependence of the catalytic 
bodies upon the NS species. Values of these constants for the 11-species air model are 
given in table I. 

The reaction rates in equation (1) or equations (2c) and (2d) are expressed in the 
Arrhenius form as: 


where 


k f , r = A f r T Bf,r exp J T ) , 
k b . r = A bir T Bb 'exp(-T D JT ) , 


s 

± 

s 


mole 

3 

cm 

mole 

3 

cm 



(3a) 


(3b) 


NJ 

ctr = E a 1>r - 1 

i=l 

NJ 

Pr= IPi.r-1 

i=l 


(4a) 

(4b) 


and T 0fr and T Dbr are the characteristic reaction temperatures for the forward and 
backward reactions, respectively. Values for the reaction rates kf r and k b r are tabu- 
lated in table II for the 11-species air model. For a specified temperature, density, and 
species composition, equations (2) through (4) can be used to obtain the production 
rate of a species i in a multicomponent gas by employing the catalytic body efficiencies 
and reaction rates from tables I and II. The first seven reactions and reaction rates in 
table II are taken from reference 16 and were employed in reference 13 for the 7- 
species air model (N ,0 ,N Z ,0 2 <N0 ,NO + ,e~). Reaction rates for reactions 8 through 
20 of table II are taken from reference 17. Some of these reactions have been 
regrouped here (and in ref. 18) through the use of third bodies M t through M 4 which 
is similar to the approach in reference 13. 

The reaction rates given in table II were originally used by Blottner (ref. 13) and 
Dunn and Kang (ref. 17) in the context of a single temperature assuming thermal 
equilibrium. One may use Park’s (ref. 19) guidelines (in the context of his two- 
temperature model) for defining the rate controlling temperature in dissociation and 
electron impact ionization reactions under thermal nonequilibrium conditions. This 
was done in reference 20, for example. 
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Park has suggested the use of a temperature weighted with the vibrational tem- 
perature to characterize dissociation reactions, based on the preferential dissociation 
concept. The reaction rates in Park’s model are assumed to be dictated by the 
geometric average temperature 

T av — s/ T^Ty,^ (5) 

and the dissociation reaction rates are given by 

*/.r = Af.JaJ' CX P(- T D,'JTav) (6) 

Treanor and Marrone (refs. 21 and 22) have suggested a more rational (but slightly 
more difficult) way than Park’s to account for the effect of vibrational relaxation on 
dissociative reactions with the preferential dissociation concept. They have suggested 
the use of a vibrational coupling factor (refs. 20,21, and 22) with the dissociation reac- 
tion rates obtained under the assumption of thermal equilibrium. 

Recent work of Jaffe (ref. 23), based on collision theory and using methods of 
statistical mechanics, found no evidence of preferential weighting to any particular 
energy mode in obtaining the total energy available in a collision whether it is an elas- 
tic, inelastic, or reactive encounter. Jaffe found that the multi-temperature effects on 
the reaction rates were small for dissociation. These findings were supported by those 
of Moss et al. (ref. 24), who carried out flowfield analyses employing the Direct Simu- 
lation Monte Carlo (DSMC) approach. Thus, a weaker dependence of k f <r on T vib 
(such as the one recently suggested by Sharma et al. (ref. 25) as T av = T^T^) might 
be more realistic, especially for highly energetic flows. 

It is obvious that the multi-temperature kinetic models for high energy flows 
based on both preferential and non-preferential dissociation assumptions employ some 
degree of empiricism. They exemplify the degree of uncertainty which exists in 
modelling the multi-temperature kinetics. Quantum mechanical studies of the type in 
reference 23 supplemented by non-obtrusive laser diagnostic studies would be desir- 
able to establish these models on a sounder basis. 


Species Thermodynamic Properties and Mixture Formulas 

Thermodynamic properties (i.e., C pi and h t ) are required for each species con- 
sidered in a finite-rate flowfield calculation. For calculations with chemical equili- 
brium, the free energies F t are also required. Since the multicomponent gas mixtures 
are considered to be mixtures of thermally perfect gases, the thermodynamic proper- 
ties for each species are calculated by using the local static temperature. Then, proper- 
ties for the gas mixture are determined in terms of the individual species properties 
through the relations: 

NS 

A=L C.h, (7a) 

*= l 


with 


T 


hi — J C p [dT 

298 


(A h/) Tm 


298 K 


(7b) 


and 



(8a) 
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NS 

Cp, ~ 2 CiCp.i 

l-l 


with 


r ^ 



(8b) 


Note that C Pf for a mixture defined through equation (8a) is the ’frozen’ specific 

heat. This definition does not account for species production or conversion due to 
chemical reactions. Frozen specific heat is commonly employed in defining the Prandtl 
and Lewis numbers! for a mixture and is related to the mixture enthalpy h through 
the relation 

NS 

= 1,0^ (9) 

p i-i 



Expressions for C p i using the partition function approach were obtained in refer- 
ence 26, whereas references 27 and 28 employed a virial coefficient method. The par- 
tition function formulation is quite accurate at low temperatures; however, this formu- 
lation becomes less accurate at high temperatures because of the introduction of errors 
from several sources such as non-rigidity of the rotor and anharmonicity of the oscilla- 
tor. The virial formulation defines the potential-energy surface more accurately at 
higher temperatures in comparison with the partition function approach. Using the 
virial formulation, Browne (refs. 27 and 28) obtained the thermodynamic properties as 
corrections to those of the monatomic gas in terms of the first and second virial 
coefficients and their temperature derivatives. These data can be used to separate the 
contributions of different internal energy modes to the specific heat as is required in 
the case of multi-temperature flowfield models. For example, in a two-temperature 
model (refs. 19 and 20), advantage may be taken of the fact that the translational and 
rotational energy modes are fully excited at room temperature, and therefore, the heat 
capacities for these modes are independent of temperature. The combined 
vibrational-electronic specific heat for species i, (C Pti ) ve , can then be evaluated by 
using the value for the total specific heat, C p i , evaluated at temperature T ve and sub- 
tracting out the constant contribution from the translational and rotational specific 
heats. This can be described (see fig. 2 for the various contributions) by 

l(C pJ ) ve ) T " = [C Pti ] Tm - (Cp,i) lr - (C pJ ) rol (10) 


The enthalpy, h it for a two-temperature model can be evaluated similarly, since 
contributions from the translational and rotational modes are linear with temperature. 
Therefore, the vibrational-electronic enthalpy for species i, (/t ( ) ve , can be obtained 
from the specific enthalpy h, evaluated at temperature T ve and subtracting out the con- 
tribution from the translational and rotational enthalpies evaluated at T ve as well as the 
enthalpy of formation. 

[(^)vJr w = [hi] Ty . ~ [<.C p ,i)„ + ( C pJ ) rot KT ve - T rtf ) - (A h/) Tnf (11) 


t See Appendix A for the various definitions. 
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The specific enthalpy from all the contributions of internal energy modes can 
then be obtained by adding the contributions of translational and rotational enthalpies 
(evaluated at the translational-rotational temperature) and the enthalpy of formation 
to the vibrational-electronic enthalpy as follows 

hi(T,T w ) = [(A,)* ] Tyt + [(C pJ ) lr +(C pi ) rol KT-T ref ) + (A h/ ) Tnf (12) 


Reference 26 has provided curve fits for C p i and h t , whereas references 27 and 
28 have provided tabulations of these data. Since the use of curve fits reduces the 
expense of computing the original functional relations, these thermodynamic proper- 
ties have been curve fitted here as a function of temperature for the temperature 
range of 300 K < T < 30000 K. These curve fits include those of reference 26 for 
the range 300 K < T < 6000 K, and new curve fits are provided to the tabulated 
values of references 27 and 28 for the range 6000 K <, T < 30000 K. The following 
polynomial equations are employed for these curve fits: 

Specific heat: 

q . 

— = A x + A 2 T + A jT 2 + A 4 T 3 + A 5 T 4 (13) 

^ univ 


Specific enthalpy: 


A« 

R T 

lx univ 1 


~ A i + 


a 2 t 

2 


+ 


a 3 t 2 a 4 t 3 a 5 t 4 a 6 

— - — + + — - — + — 

3 4 5 7 


(14) 


For equilibrium calculations, the following curve fit for the free energies F, may 
be used: 


R T 

lx univ 1 


= A x [\-ln{T)) - 


A ?T A -iT 1 


A 4 T 3 

IT 


A£*_ A 6 - Aj 

20 + T 


(15) 


where F°‘ is the free energy of species / at 1 atm pressure (standard state). 

A tabulation of the polynomial constants (A t to A 7 ) for the 11-species air model 
is given in table III. A least squares technique has been used for obtaining the poly- 
nomial coefficients for the five temperature ranges between 300 and 30000 K. To 
assure a smooth variation of thermodynamic properties over the entire temperature 
range, values of A ! to A 7 should be linearly averaged across the curve-fit boundaries 
(i.e., 800< T < 1200, 5500< T < 6500, 14500< T < 15500, and 24500< T < 25500). 
An example subroutine which evaluates the polynomial curve fits and performs the 
linear averaging is presented in appendix B. This routine may be easily modified to suit 
the user’s requirements. 

It may be mentioned here that temperature, T, in equations (13) through (15) is 
in °K . With the universal gas constant in call gm- mole- K , the specific heat and 
enthalpies will have the units of call gm- mole- K and call gm- mole respectively. The 
constant A 6 in these equations is related to the heat of formation through the relation 

A 6 R univ = (AA/)r=298tf (16) 


Figure 3 compares the values of specific heat obtained from the polynomial curve 
fit (eq. (13)) with the data of McBride et al. (ref. 26) and Browne (ref. 28). Values 
provided by Hansen (ref. 29) are also included for comparison. Hansen’s values begin 
to deviate from those of references 26 and 28 beyond 4000 K for 0 2 (fig. 3(a)), 
beyond 8000 K for N z (fig. 3(b)), and beyond 6000 K for NO (fig. 3(c)). This may 
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be due to the rigid-rotator and harmonic -oscillator partition function employed by 
Hansen. At higher temperatures, the population of levels corresponding to the non- 
parabolic regions of the potential energy curve are no longer negligible, and it becomes 
necessary to introduce the non-rigidity and anharmonicity corrections into the energy 
levels of the molecules as was done in reference 28. 


Species Transport Properties and Mixture Formulas 

The transport properties required in flowfield calculations are viscosity, thermal 
conductivity, and diffusion coefficients. The collision cross sections required for these 
properties have been recomputed herein using the same molecular data used previ- 
ously by Yos (refs. 7,8, and 9). The computational techniques employed in the calcu- 
lations are described in references 30 and 31 which give details of the NATA (None- 
quilibrium _Arc Tunnel Analysis) code (ref. 32). In NATA, the average collision cross 
sections for the collisions between species i and j are calculated from basic 

cross section data as functions of temperature and gas composition for each pair of 
species in the mixture. The basic data are either in tabular form or are given as simple 
analytical functions of temperature or composition. NATA contains twelve methods 
or options for calculating the cross sections re Cl^ 2 ' 2 \ and B * (the ratio of 

cross sections). The options include using the Coulomb cross section for the electrons 
and ions plus exponential potential (ref. 33) and Lennard-Jones (6-12) potential (ref. 
34) for neutral species in high and low temperature ranges, respectively. The formulas 
employed in NATA to compute the transport properties from the collision cross sec- 
tions are obtained from an approximation (ref. 35) to the rigorous first-order 
Chapman-Cowling expressions. 


Single Species Transport Properties 


The viscosity jx t and frozen! thermal conductivity i of a gas containing a sin- 
gle molecular species are given, to a good approximation, by the formulas (ref. 34, 
Chapters VII and VIII): 


5 V KiriikT 

16 iznW ’ 


$rn 

cm - sec 


(17a) 


= 2.6693x 10" 5 


rc V MjT 


cm - sec 


(17b) 


with 


K f , = K lr>i + K 


k mt y i 


K lr ,i = 


75 W k / cT I mi 

64 rt nr- 2) ’ 


erg 

cm - sec -K 


= 1.9891X 10 -4 


~^F r ’ 


cal_ 

cm - sec - K 


15 Mum v 
4 


cal_ 

cm-sec-K 


(18a) 

(18b) 

(18c) 

(18d) 


t See Appendix A for various definitions 
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and 


K = — 
int ' 1 8 




k^ nkT I mi 




,_ s rcVr/M, 


'p.‘ 


cm- sec- K 

cal 


= /?, 


M-i 

A*i 


*m u> 

^ ^ univ 


2 J 

cm-sec-K 

' P D u ' 


7,i 

_5 


cal 

P, J 

[ R univ 

2 

» 

cm-sec-K 


(18e) 


(180 


(18g) 


where A’, r<( is the translational component of the thermal conductivity and K int%i is the 
component of thermal conductivity resulting from the diffusion of internal excitation 
energy of the molecules. In equations (17a) through (18g), is the molecular mass 
in gm, k is Boltzmann’s constant (1.38x 10“ 16 erg/K), T is the gas temperature in K, 
rcQ’( /,/) are average collision cross-sections for the molecules in A 2 , (where 
1 A = 10 _8 cm), Mi is the molecular weight in gmlgm-mole, 
Runiv- 1.987 call gm- mole- K , C pi is the specific heat at constant pressure in 
call gm- mole- K , p is the density in gm/cm 3 , and D u is the coefficient of self 
diffusion in cm 2 / sec. 

From equations (18a), (18d), and (18g), the frozen thermal conductivity may 
also be expressed as 


*/, = 


*««ivHi 
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(19) 


The factor (pD a /|x { ) appearing in equations (18g) and_(19) is the reciprocal Schmidt 
number and is related to the collision cross sections nTTtf'D through the relation (ref. 
34) 

pp, 6 

h, 5 Km 1 ' 1 * 


( 20 ) 


with D u defined as 


Du = 


= 2.6280x 10" 3 


3 V KmikT 

I k£T( u > 

W T^/Mi 


cm 

sec 


*m u) P 


cm^ 

sec 


(21a) 


(21b) 


where p is the pressure in atm . The coefficient of self-diffusion D u , must be regarded 
as somewhat artificial. It is more correct to regard it as a limiting form of the 
coefficient of binary diffusion. The ratio ( 2 ’ 2) /rc£2 appearing in equation (20) 
is a very slowly varying function of temperature, T , and hence p is very nearly 
constant. This factor appearing in equation (19) has a value of unity in the Eucken 
correction (ref. 34) and is close to 1.32 for the Lennard-Jones potential over a large 
temperature range (ref. 36). 

The collision cross sections or , in general, are the weighted 


t The collision cross sections designated here as 
given in reference 34 and as £2 given in reference 31 


are the same as 


7tc 2 Q('' 


<y 
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averages of the cross sections for collisions between species i and j. These have been 
defined (refs. 7,8 and 9) as 


> = 


oo n 

j j exp(-y 2 )y 2s +' i (l-cos l x)4KCS i jSinxdxdy 
0 0 

oo Jt 

J J exp(-Y 2 )Y 2 j '*‘ 3 (l-cos'x)sinxdx^Y 

o o 


( 22 ) 


where a,y = CTy(x,g) is the differential scattering cross section for the pair (/,/), x is 
the scattering angle in the center of mass system, g is the relative velocity of the col- 
liding particles, and y = ([m,m y /2(m,+m y )Jfc7] l/2 )g is the reduced velocity. For colli- 
sions between the similar species, equation (22) yields jc TTff’ 1 '* required in equations 
(17) and (18). 

The collision cross sections n. and nTT' L [- 2,2) employed here are the same as 
those used by Yos (refs. 7,8,9). The cross sections for the neutral species 
N ,N 2 ,N0,0 , and O 2 were taken from the tabulations of Yun et al. (ref. 37), for 
temperatures up to 15,000 K . Above 15,000 K , the cross sections for atomic N were 
obtained by extending Yun’s calculations to 30,000 K using the same input data and 
techniques as were used in his work. The cross sections for the remaining species 
N 2 ,NO ,0 and 0 2 were extrapolated to roughly 30,000 K assuming the same tem- 
perature dependence as calculated for N . 

For the ionized species, the calculations used effective Coulomb cross sections 
chosen to make the computed transport properties for a fully ionized gas agree with 
the correct theoretical results (ref. 38) as discussed in references 7 and 9. The specific 
formulas used in the calculations are 


n a ( J' 2) - 1.29<2 C for electrons (23a) 

7crr// 2 ’ 2) = 1.36Z 4 (2 c for ions (23b) 

j * 1 ' 13 = tc Q = 0.795 Z 4 <2 c for electrons and ions (23c) 


where Z = 1 for singly ionized species and 


Qc = 

The shielding parameter A is defined as 


(IcT) 2 


In A 
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2.09x 10"‘ 
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+ 1.52 


2/3 ' 


1/2 


io n Pe 


(23d) 


(23e) 

(230 


where T is the temperature in K and p e =n t kT is the electron pressure in atm. In 
equation (22d), e=4.8x 10 _10 es« is the electron charge. Equation (23) is applicable 
for electron pressures of the order 1 atm or below and becomes less accurate as the 
electron pressure increases above atmospheric (ref. 38). 

In the present report, the transport properties of the ionic species are provided 
for a nominal electron pressure of p e = 1 atm. For any other electron pressure, it will 
be necessary to correct the tabulated transport properties of the ionic species according 
to the formula 
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MjV) _ K f,i(Pe) _ ^jrAPe) _ ^mt/Pe) _ In A( 1 atm) 
11,(1 atm) Kf j(latm) AT lrl (l atm) £*,,(! atm) InA (p e ) 
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(24a) 


Similarly, the collision cross sections nCTf} ,s) (for the pair of species where both are 
ions or electrons or a combination of the two) for any other electron pressure p e may 
be obtained from the values provided herein for p e = 1 atm by employing the relation 
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Equation (24b) is also applicable for i= j (single species). 

In calculating the contribution K inti of the internal energy states to the thermal 
conductivity for the atomic species N and O, the diffusion cross section in 

equations (18e) or (18f) has been set equal to the corresponding charge-exchange 
cross section for the atom and atomic ion (refs. 8 and 39). As discussed in reference 
7, this approximation allows for the effects of excitation exchange in reducing the con- 
tribution of internal energy states to the thermal conductivity in a gas of identical 
atoms. 


The individual species viscosities and thermal conductivities computed using 
equations (17) and (18) have been curve fitted herein as a function of temperature by 
employing the following relations 

jx,- = * 'T ‘ * ** ' , g” (25) 

cm- sec 


K /, = 


Eg [*K V*T) 3 +B K ( lnT) 2 + 
e K f.iJ K fS K f.i 


C K In 7* 

K /,i 




cal 

cm- sec -K 


(26) 


The curve-fit coefficients appearing in equations (25) and (26) are given in tables 
IV and V for the 11 chemical species. These coefficients yield values of the viscosity 
and frozen thermal conductivity of ionic species at a nominal electron pressure of 
1 atm and should be corrected for any other electron pressure (or number density) by 
using equation (24a). 

Figure 4 displays typical results from the viscosity curve fit (eq. (25)) of equation 
(17) for some of the neutral and charged species. The frozen thermal conductivity 
curve fit (eq. (26)) of equations (19) through (21) is shown in figure 5 for the same 
species. It may be noticed from these figures (and in tables IV and V) that a higher 
order curve is needed for the charged species because of the inflection of the viscosity 
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and thermal conductivity curves at low temperatures. Higher order curve fits are also 
needed for the thermal conductivity of the neutral molecular species. 


Transport Properties of Multicomponent Mixtures 

Rigorous kinetic theory formulas which have been derived directly from a solu- 
tion of the Boltzmann equation using the classical Chapman-Enskog procedure (refs. 
33 and 40) are available for obtaining the transport properties of a gas mixture from 
the molecular constituent species. In the first Chapman-Enskog approximation, for- 
mulas for both the viscosity and translational component of thermal conductivity, 
K lri , of a gas mixtures are of the general form 


n 


lv 


1 vl • 


[H (1) or 


*1 


mixture 


x v I 0 


\A tJ \ 


(27) 


where x- t is the mole fraction of the ith species, v (2 NS) is the total number of 
species present in the mixture, and the matrix elements can be expressed in the 
form 

V 

Ajj Aji x | x j a ( y + hij{XiAi + x {X 1 Qj ) (28) 

/= 1 

where a,y=ay,- and 5,y is the Kronecker delta. Elements A, and a,y are defined subse- 
quently. The superscript 1 on p. or K tr indicates that equation (27) is the first 
Chapman-Enskog approximation for the transport property. Further details for obtain- 
ing the transport properties by employing the first Chapman-Enskog approximation are 
given in reference 35. In principle, the problem of calculating the transport 
coefficients for a given mixture consists of two parts: first, the determination of the 
collision cross sections JtQ, y for all possible pairs of species (/,/); and second, the 
evaluation of the Chapman-Enskog formulas. The amount of computation required to 
evaluate the mixture transport properties is greatly reduced if approximations to the 
complete Chapman-Enskog formulas are employed. Reference 7 has provided approx- 
imate formulas for the transport properties based on the relations developed in refer- 
ences 41, 42, and 43. Earlier, references 10, 11, and 44, and more recently, refer- 
ences 12 and 45 have also provided approximations to the Chapman-Enskog formulas. 
These are apparently the most satisfactory of the many simplifying approximations for 
the mixture viscosity and thermal conductivity which have been suggested by various 
authors. However, effects of the elements A^ in the Chapman-Enskog formula are 
completely neglected in Brokaw’s approximation (ref. 42) so that this approximation 
always gives too large a value for the transport properties. In the Buddenberg-Wilke 
(refs. 10 and 44) and Mason-Saxena (ref. 11) formulas, the effects of these elements 
are accounted for by means of a single empirical constant which is assumed to be the 
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same for all gas mixtures. The approximation used by Peng and Pindroh (ref. 41) 
represents an attempt to take account of the nondiagonal elements explicitly to first 
order at the expense of a somewhat increased calculational effort. Armaly and Sutton 
(refs. 12 and 45) neglected the nondiagonal matrix elements A ,y in a manner similar to 
Brokaw (ref. 42). However, they did not force the value of A-j, defined as 


A ij 


a 


( 2 . 2 ) 


0 (U) 


(29) 


to be equal to 5/3 and 5/2 in their approximations for viscosity and thermal conduc- 
tivity, respectively.! They assigned different values to A* for ion-atom and neutral 
atom-molecule interactions. From the computer time, storage, and simplicity point of 
view, references 10 and 11 would appear to be adequate for non-ionized gas mixtures, 
whereas references 12 and 45 would be useful for computing the viscosity and transla- 
tional component of thermal conductivity for an ionized gas mixture. 

In all the approximations to the Chapman-Enskog formulas for viscosity and 
translational thermal conductivity discussed thus far, the transfer of momentum or 
energy from one species to another by collisions has been either neglected or has been 
accounted for by an empirical constant. This transfer process, which is represented by 
the nondiagonal elements A, ; - in the Chapman-Enskog formula (eq. (27)), has the 
effect of making the less conductive species in the mixture carry a larger fraction of 
the transport. This process, therefore, reduces the overall conductivity of the mixture 
below the value which it would have if the transfer process were neglected. 

In reference 35, Yos obtained approximations to the Chapman-Enskog formulas 
which account for the effects of the above mentioned transfer process between 
different species. These cross sections reproduce the results of the first Chapman- 
Enskog formula (eqs. (27) and (28)) to within a fraction of a percent for all cases con- 
sidered in reference 35 and are also simpler to use than the latter. Based on the cross 
sections developed by Yos (ref. 35), the following formulas may be used to compute 
viscosity and the translational component of thermal conductivity 

NS 

%Xi/(Ai + a av ) 

[H (1) or K?'] = — (3°) 

1 - a av '£x i l(A i + a av ) 

i=i 


Here, NS is the number of species in the gas, x t is the mole fraction of the t'th 
species, and a av is an average value of the nondiagonal matrix elements a {j contained 
in equation (28) which is defined as 
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I XiXj 
• J m 1 
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NS 
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(31a) 


with 


t With A* = 5/3 and 5/2, A tJ becomes identically equal to zero in the Chapman-Enskog 
formula. 
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NS 

Ai = iMi 

/-i 


(31b) 


For the viscosity, the quantities a and A { j appearing in equations (31a) and (31b) are 
defined as 

a - = — r ?A 01— 

Mi 

For translational thermal conductivity, the above quantities are defined as 


(32a) 

(32b) 
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(33a) 


(33b) 


In these equations, M, is the molecular weight of the ith species, 
N A = 6.0225x 10 23 molecules/ gm- mole is Avogadro’s number, and 

* = 1.3805x 10 -16 erg/K is Boltzman’s constant. The remaining quantities are defined 
as 
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(34) 

(35) 

(36) 


and the collision integrals are weighted averages of the cross sections defined 

through equation (22). 

The frozen thermal conductivity K f employed in defining the frozen Prandtl and 
Lewis numbers! can be obtained from the modified Eucken approximation (ref. 43) 

(37) 


K f = K lr + K h 


tnt 


where K tr is the translational component of the thermal conductivity given by equa- 
tions (30) to (36) and K inl is the component of thermal conductivity resulting from 
the internal excitation energy of the molecules, defined (ref. 7) as 


NS 

K inl = 

i=l 
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kN A 2 


Xi 


NS 


( 38 ) 


+ See Appendix A for various definitions 
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Here, C pi is the specific heat at constant pressure of the ith species. 

The binary diffusion coefficient D t y needed to obtain the binary and multicom- 
ponent Lewis numbers (refs. 13, 14, and 34) is obtained from the complete first 
Chapman-Enskog approximation (ref. 34): 

Du = — - = ttt , cm 2 1 sec (39a) 

P 

with 




kT 

A ,0) 

“y 


cm -atm! sec 


(39b) 


Here p is the pressure in atmospheres. Equation (39b) for D,-, has been curve fitted 
in the present work by using the expression 


_ D 5 [A p . (In T y 2 + B B In T + CjJ 
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(40a) 


where T is the temperature in K. The curve-fit coefficients of equation (40a) are 
given in table VI for the different interactions occurring in the 11-species air model. 


The values of Z> ( y obtained from equation (40a) with the curve-fit coefficients 
given in table VI are for a nominal electron pressure of 1 atm. If the pair of interact- 
ing species are both ions, both electrons, or any combination of the two, then D it may 
be corrected for any other electron pressure by employing: 
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(40b) 


If the approximations of references 41 through 43 are employed in place of the more 
exact formulas given by equations (30) through (36) t , then the following approxi- 
mate formulas for the mixture viscosity and thermal conductivity may be used (ref. 
7): 
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where A- 2) is given by equation (35) and a,y is defined as 


a v = 1 + 


[1- (M,/Af ; )][0.45- 2.54 (Af.VAfy)] 

[1 + (Af,/Af ; )] 2 


(41) 


(42a) 


(42b) 


+ The approximations of references 41 through 43 may be valid only for the range of condi- 
tions for which they have been developed and not for general application because of the very ap- 
proximate analysis of the nondiagonal matrix elements Ay of equation (27). 


- 19 - 


To obtain the viscosity n and the translational component of thermal conductivity 
K tr of a gas mixture from either the more exact equations (30) through (36) or from 
the approximate relations in equations (41) through (42), the binary collision cross 
sections kTi^ , s) and their ratios are needed in equations (34) through (36). These 
cross sections, defined by equation (22), are the same as those used by Yos (refs. 7,8, 
and 9). In addition to the formulas given by equation (23) for the Coulomb collision 
integrals, the following relations are also employed: 

izQ ( J a) = 0.795Z 4 Q c (43a) 

KCT e ( , 2 ' 2) = 1.29(2, (43b) 

where the subscripts e and I represent electrons and ions, respectively, and Z and Q c 
are defined in equation (23d). 

The NATA code (refs. 30,31, and 32) employed to obtain the collision cross sec- 
tions, KCl[j l ' s \ contains default provisions for estimating some cross sections if they 
are not specified explicitly in the built-in database or the input. The defaults are sum- 
marized as 

(i) If both species are ions, then the Coulomb cross sections given by equations (23) 
and (43) are used. 

(ii) If one species is neutral and the other ionized, then the formula 

%SZlj’ s) = a (Ls) T~ oa (44a) 

is employed with the constants defined in the code. 

(iii) If both species are neutral and are unlike, then the cross sections are estimated 
using the simple mixing rule 

KCTfj’ s) = j(yj it 4 V Tclfp 7 ) 2 (44b) 

The built-in data in NATA specify steps for calculating the cross sections for the 
like-like interactions of ten species (e~ ,N 2 , 0 2 ,N , O ,NO , NO + , N + , 0 + , andN £ ) and 
for those unlike interactions for which experimental or theoretical cross sections are 
available in the literature. The cross sections for O £ are taken to be the same as 
those for in these calculations. NATA contains twelve methods or options for 
calculating 7tfl^ 1,1 \ k £2 ( y 2,2 \ and B*j which are described in detail in references 30 
through 32. Note that the accuracy of the calculated transport properties is largely 
determined by the accuracy of the input data for the cross-section integrals. A brief 
discussion of the values employed in the present study is given in reference 7. Refer- 
ence 46 provides a comparison between theory and experiment for the thermal con- 
ductivity of nitrogen up to temperatures of 14000 K . Fairly good agreement between 
the two is shown there. 

Figures 6 and 7 give viscosity and frozen thermal conductivity values, respec- 
tively, obtained for equilibrium air at 1 atm by employing the presently computed col- 
lision cross sections for the constituent species. The present calculations employ equa- 
tion (41) for the calculation of mixture viscosity and equations (37), (38), and (42) 
for the mixture frozen thermal conductivity. In addition, the collision cross sections 
used are nearly identical to those obtained in reference 9. The mixing laws employed 
for the viscosity and thermal conductivity shown in figures 6 and 7 are accurate for 
sub-ionization temperatures (less than 9000 K at 1 atm) only. For ionized flows, more 
accurate mixing laws of the type given by equation (30) should be employed. Figures 
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6 and 7 also include viscosity and thermal conductivity values from other researchers 
for comparison. In these figures, the predictions of Peng and Pindroh (ref. 41) and 
Esch et al. (ref. 47) are based on the Buddenberg- Wilke (ref. 44) type of mixture law 
and have the same level of approximation as the present calculations. The mixture 
laws employed by Hansen (ref. 29) contain a somewhat lower level of approximation. 
In Hansen’s work, the viscosity is computed by using the simple summation formula 
for a mixture of hard sphere molecules, whereas a linearized expression with Eucken’s 
assumption (ref. 34) is employed for the frozen thermal conductivity. The viscosity 
values obtained in the present work are in good agreement (fig. 6) with the values 
obtained by Peng and Pindroh (ref. 41) and Esch et al. (ref. 47), presumably due to 
the similar mixing laws employed. Hansen’s (ref. 29) predictions are lower, especially 
at higher temperatures. Even with the beginning of dissociation of molecular nitrogen 
at about 4000 K (when the dissociation of molecular oxygen is almost complete), 
Hansen’s values of viscosity are not much different from the Sutherland’s law. For 
the frozen thermal conductivity (fig. 7), the present values are in agreement with 
those obtained by Peng and Pindroh (ref. 41) up to temperatures of about 9000 K. 
The values obtained by Esch et al. (ref. 47) deviate from the present values beyond 
6000 K. This may be due to the constant cross-sectional values employed for the ion- 
ized species in the 8000 to 15000 K temperature range in reference 47. Hansen’s 
predictions of thermal conductivity are lower than the other data and are closer to the 
Sutherland values up to temperatures of about 4000 K. Again, the differences 
between the present computations and those of Hansen are presumably due to the 
somewhat more rigorous mixing laws employed herein. 

There are 121 possible binary interactions for the dissociating air having eleven 
species. Therefore, 121 values of each of the collision cross sections kTT^- 2) 

and the collision cross-section ratio B*j are required to evaluate the transport proper- 
ties. By use of the symmetrical equality, i.e. (i,j) = O',/), only 66 values of each cross 
section and cross-section ratio are required. These values have been curve fit in the 
present study as a function of temperature for a nominal electron pressure of 1 atm 
using the following relations 


jcQ< U) =e * T ■' " # , 

A 2 

(45) 

KClf' 2) =e ■' T ■' ■' *' , 

A 2 

(46) 

D . [A .(InT ) 1 + B . + C ,] 
B*j = e "T “ “ “ 


(47) 


For electron pressures different from unity, the formula given by equation (24b) 
is used to correct the cross sections for the ionic species. No such correction is 
required for the cross-section-ratio parameter B*j. 

Curve-fit coefficients appearing in equations (45) through (47) are given in Table 
VII through IX for the various interactions in an 11-species air model. Figure 8 illus- 
trates some typical curve fits obtained by employing equation (40a) and equations (45) 
through (47) with the associated constants. The figure compares the computed values 
of binary diffusion coefficient, collision integrals, and collision integral ratio with the 
resulting curve fit for different interaction pairs of neutral and ionized species, includ- 
ing electrons. The collision integral ratio, B*. is almost constant with temperature as 
shown in the figure and was fitted with the lower order curve fit where possible. How- 
ever, the collision integrals TtQ-}' s) for the charged species pairs exhibit more complex 
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behavior and require higher order curve fits. 

The transport properties of a multi-temperature gas mixture may be obtained by 
following the approaches of references 20 and 48. These references have used equa- 
tions (41) and (42) for the calculation of mixture viscosity and translational thermal 
conductivity. The collision integrals for heavy particles in these equations are 
evaluated by using the heavy particle translational temperature, whereas those for elec- 
trons with any other partner are obtained by using the electron temperature. With 
these modifications, equations (41) and (42) can be written for an 11-species gas mix- 
ture as 
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(49) 


where a y is still obtained from equation (42b). Note that the above definition of K* r 
does not include contributions due to electron-heavy particle and electron-electron col- 
lisions. These contributions, defined by K e , are given later. The approach outlined 
here may also be used to obtain the mixture viscosity and thermal conductivity from 
the more exact expressions in equations (30) through (36). Equation (38) for the 
internal thermal conductivity needs to be modified for the multi-temperature formula- 
tion. The contributions resulting from the excitation of different internal energy 
modes can not be lumped together into a single term K mt for such a formulation. 
Further, equation (37) no longer can be used to obtain a frozen thermal conductivity 
for the mixture. The relation given in equation (37) may be used to obtain a frozen 
thermal conductivity with only the rotational mode contributing to the internal energy 
at the translational temperature. In general, there are four components of the internal 
thermal conductivity, similar to the molecular specific heat (fig. 2). Using these com- 
ponents, the &th component of the overall heat flux vector can be expressed as 


q k = -( K ; r + 

OX 


dT v 

dx k 



dT NS 

f + Z 9iWi 
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(50a) 


where K* r is the translational thermal conductivity defined previously by equation (49) 
and K rot , K vib and K e[ are the rotational, vibrational, and electronic thermal conduc- 
tivities, respectively, associated with these internal energy modes (ref. 49). Also, K e 
is the thermal conductivity of electrons, T is the translational-rotational temperature, 
T v is the vibrational temperature, T el is the electronic excitation temperature, T e is the 
electron temperature, x k is the &th component in a general orthogonal coordinate sys- 
tem, and the last term is the diffusive component of the heat flux vector. Further 
details of the heat flux vector q k and other definitions are given in appendix A. Two 
other heat flux vectors are required in describing the conservation equations for the 
electron energy and the vibrational energy (ref. 48). The electron heat-flux vector q k 
and the vibrational heat-flux vector q k are defined in appendix C. The thermal con- 
ductivities required to evaluate these flux vectors are obtained from the electron- 
electron and molecule-molecule collisions. These thermal conductivities are also 
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defined in appendix C. In the diffusive heat flux component, p,-, /t,, and V k 
density, enthalpy and diffusion velocity of species i, respectively. For 
temperature model, equation (50a) may be written (refs. 19 and 20) as 
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Different components of the internal thermal conductivity appearing in equation 
(50) can be evaluated from equation (38) by appropriate modifications for these com- 
ponents. For example, K rol can be obtained from 
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for partial excitation of the rotational internal energy mode if the temperature is less 
than that needed to excite the vibrational energy mode (fig. 2). Values of specific heat 
at constant pressure, C pi , appearing in equation (51a) can be obtained from the 
curve-fit relation of equation (13) by employing the rotational temperature, T rol , if 
different from the translational temperature T. When the vibrational mode begins to 
excite, the rotational mode is fully excited and equation (51a) becomes 
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2>;A tj l \T rol ) 
7-1 


(51b) 


Similarly, expressions for K vib with partial and full excitations of the vibrational energy 
mode may be written as 


(K vlb ) Partial 
Excitation 




Mi C p i (T vib ) 7 


NS 

i-l 


N a k 2 

k. J 



NS 

ZXjA[/\T vib ) 



(52a) 


and 


NS 

b) Full - 

Excitation i- 1 


Xi 


NS 

ZxjA^iT^) 

7-1 


(52b) 


Once again, the value of C p , appearing in equation (52a) can be obtained from the 
curve-fit relation of equation (13) by employing T vib . Also, the vibrational energy 
mode is fully excited when the electronic contribution becomes significant (as shown 
in fig. 2 for the specific heat at constant pressure). 

The rotational and vibrational energy modes are almost fully excited at their 
respective characteristic temperatures. Therefore, for rotational and vibrational tem- 
peratures greater than these characteristic temperatures, equations (51b) and (52b) can 
be employed to obtain K roi and K vlb , respectively. 
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as 


Expressions for K ei with the excitation of electronic energy mode can be written 


*el = 


NS 


*1 


M i Cp,i(T el ) jj 

N a k 2 

NS 

ZXjAW(T el ) 
j= 1 


(53) 


In equation (53), higher electronic degeneracy levels begin to contribute to K el (and 
C pi , etc.) as the characteristic temperatures for the electronic excitation of those 
degeneracy levels are reached. 

Finally, the thermal conductivity for free electrons, K e , appearing in equation 
(50) may be obtained from the modified form of equation (42) 


K e 


ii 

4 NS 

Z L45xjA e f(T e ) 
/= i 


(54) 


This thermal conductivity results from the collisions between electrons and other 
species, including other electrons. 


Concluding Remarks 

The present work provides the reaction rate coefficients and thermodynamic and 
transport properties for an 11-species air model which are needed in analyzing the high 
energy flow environment of currently proposed and future hypersonic vehicles. 
Approximate and more exact formulas have been provided for computing the proper- 
ties of partially ionized gas mixtures in chemical and thermal nonequilibrium around 
such vehicles. 

The work presented here uses the best estimates of available data needed to com- 
pute properties of the 11-species air model. However, there is need for improved 
data, especially for air in thermochemical nonequilibrium. There is also some degree 
of uncertainty about reaction rate coefficients at high temperatures. Both the theoreti- 
cal and experimental bases of the multi-temperature kinetic models need to be 
strengthened. The virial coefficient approach as compared to the partition function 
analysis for obtaining the thermodynamic properties at higher temperatures appears 
sufficient but needs further verification. For transport properties, the input data for 
obtaining the collision cross sections at high temperatures are not adequate. 
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Appendix A 


Heat Flux, Frozen and Reactive Prandtl and Lewis Numbers, and Associated 
Definitions 

The *th component of overall heat flux vector q k for the dissociating and ioniz- 
ing 11-species air model in a multi-temperature formulation can be expressed asf 

. . dT ar v ST,, ST, It 

q k = -( K ; r+ K rol )^- K vU> j^- K tl ^~ (Al) 


Sx k 


Sx l 


where the various symbols are explained in the main text after equation (50a). For 
thermal equilibrium conditions, equation (Al) may be written as 


q k = -(*;+ K r „ + K ^ + K. l + K t )-p T +Zf> i h i V? 

OX I 


or, 


with 


and 


ar 11 

OX i»l 


K f = (C+^.)+^ta, 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 


K ini ~ K rot "I" ^ vib ^ tl 

It may also be noted that under thermal equilibrium conditions 

K, r = K^+K. 

where K„. is given by equation (42), K’ lr is obtained from equation (49) and K, is 
evaluated from equation (54). The last term in equations (Al), (A2) and (A3) 
represents the diffusion contribution to the heat flux vector. The diffusion mass flux 
of species i, J k , is related to the diffusion velocity V k through the relation 

■/* = PiV? (Al) 

and for the case of binary diffusion may be expressed as (ref. 50) 

* dQ 

J! pDi 'd?‘ (A8) 

when the species i diffuses into any other species j. Here, C, is the mass fraction of 
species i and D„ is the binary diffusion coefficient. Using relation (Al) and (A8), 
equation (A3) can be written as 


ST 


ii 
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3C; 

Sx k 


or, 


l k = ~K, 


Sx k 
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(A 9) 


(A 10) 


t Here we have used the sign convention that the heat flux vector is positive in the normal 
direction away from the surface. 
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or, 



where K is the total effective thermal conductivity defined as 


K = Kf + ^Dijh; 


i=i 


dQ 

ar 


(All) 


(A12) 


or, 


K = K f + K d 

with the diffusion contribution to the thermal conductivity, K d , defined as 


11 dCi 

K d — p^Djjhi 

i=i ai 


(A13) 


(A14) 


An alternate definition of the total thermal conductivity, K , is given in reference 7 as 

K = K f + K r (A15) 

where K r is the reaction thermal conductivity! defined as 


NIK 

K r = *1 
(=1 


(Aht/R^T)* 


NS NS 

I - <*,•.,)/*< H K0i., -*.,)*>- (P,./ ~ ajMAl? 

i= l /= l 


(A16) 


where k is the Boltzmann constant, NIR is the total number of independent reactions 


NS 


in the system, NS is the total number of species in the system, A h,= £(P M - a, is 


<=i 


the heat of reaction per gm-mole for the lih reaction, is the mole fraction of species 
i. A// 5 is defined through equation (34), and a,-,/ and Pi.j are the stoichiometric 
coefficients for reactants and products in the reaction given by equation (1). The ther- 
mal conductivities K d given by equation (A14) and K r given by equation (A16) are 
equivalent. 

Similar to relations (A12) and (A13), the total specific heat at constant pressure 
may be defined as 


or, from equation (9), 


or, 


C,= 


C = Cn 

p pf 


ar 


ii ac, 


C. — C- + 

p Pf Pd 


with 




» dCi 


(A 17) 


(A 18) 


(A19) 


(A20) 


+ The formula for Kr given here already includes the effects of ambipolar diffusion (ref. 51) 
on the reaction conductivity of an ionized gas. 
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where C Pj is the diffusion contribution to the specific heat at constant pressure. 

Using the various definitions for the thermal conductivity and the specific heats at 
constant pressure, the frozen and reactive Prandtl and Lewis numbers usually 
employed in dimensionless heat flux calculations may now be defined.! 

Generally, the frozen values of specific heat at constant pressure and thermal 
conductivities are employed in flowfield calculations. Associated heat flux from these 
calculations is, accordingly, expressed in terms of the frozen Prandtl and Lewis 
numbers. Alternately, one may also use total values of specific heat at constant pres- 
sure and thermal conductivity in flowfield computations. In this case, the heat flux 
may be expressed in terms of the reactive Prandtl and Lewis numbers. However, the 
total values of C p and K as well as the reactive Prandtl and Lewis numbers can be 
used only with calculations involving thermal equilibrium. This is obvious from the 
definitions of K and C p given by equations (A13) and (A19), respectively. 

The frozen Prandtl number and Lewis numbers are defined as 

Pr f = C PfV LlK f (A21) 

Le f .i = pC P/ Dij/K f (A22) 

where C P/ and K f are defined by equations (9) and (37), respectively. Generally, the 

subscript / is not used in the literature to denote the frozen Prandtl and Lewis 
numbers as is done here. 

The reactive Prandtl and Lewis numbers are defined by employing the total 
values of C p and K as 

Pr r = C p \ilK (A23) 

Le r<i = p C p DijlK (A 24) 

where C p and K are defined through equations (A18) and (A13) or (A15), respec- 
tively. Reference 29 gives definitions of the reactive Prandtl and Lewis numbers with 
K defined through use of equation (A15) and employs no subscripts on Pr or Le. 
The frozen values defined through equations (A21) and (A22) are denoted as partial 
values in reference 29 and are denoted as Pr' and Le' there. 


f Dr. Ken Sutton of NASA Langley Research Center provided these definitions in an unpub- 
lished memo. 
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Appendix B 

Sample Program to Evaluate Thermodynamic Properties From Polynomial Curve- 
Fits 

This appendix lists a sample Fortran subroutine which evaluates the species 
specific heats and enthalpies for an 1 1-species air model using polynomial curve Sts as 
functions of temperature. Five temperature ranges are used for each species for tem- 
peratures between 300 K and 30000 K. Properties evaluated near the temperature 
range boundaries are smoothed by linearly averaging the polynomial coefficients to 
assure continuous derivatives. The subroutine may be easily modified for different 
needs. 


SUBROUTINE THERMO(T.CPI.HI) 

C Computes enthalpy and specific heat for 11 species by 
C approximating polynomials. Polynomial coefficients are stored in 
C arrays A1 to A6 and are linearly averaged at the temperature range 
C boundaries. 

C 

C input: T temperature, K 

C output: CPI specific heats of the species, cal/gm-moie-K 
C HI enthalpies of the species, cal/gm-mole 

DIMENSION Al(l 1,5),A2( 1 1,5), A3( 1 1,5),A4( 1 1,5),A5( 1 1,5),A6( 1 1,5) 
DIMENSION P(6),COEF( 11,5,6) 

DIMENSION CPI(11).HI(11) 

EQUIVALENCE (Al.COEF) 

C Universal gas constant, cal/gm-mole-K 
DATA UNIR/ 1.987/ 

Coefficients are input for five temperature ranges 


K= 4 


L- 5 

IF(T.GT. 15500.)GO TO 20 
K= 3 
L= 4 


IF(T.GT.6500)G 0 TO 30 
K= 2 


L= 3 

IF(T.GT.1200.)GO TO 40 

K= 1 

L= 2 

PA= 1.0 

PB= 0.0 

EF(T.LE.300.)GO TO 50 
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PB= (l./400.)*(T-800.) 

PA= 1.0-PB 
GO TO 50 
40 CONTINUE 
PA= 1.0 
PB= 0.0 

IF(T.LE.5500.)GO TO 50 
PB= (l./l000.)*(T-5500.) 

PA= 1.0-PB 
GO TO 50 
20 CONTINUE 
PA= 0.0 
PB= 1.0 

IF(T.GE.25500.)GO TO 50 
PA= 1.0 
PB= 0.0 

IF(T.LE.24500.)GO TO 50 
PB= 0.001 *(T-24500.) 

PA= 1.0-PB 
GO TO 50 
30 CONTINUE 
PA= 0.0 
PB= 1.0 

IF(T.GE. 15500.) GO TO 50 
PA= 1.0 
PB= 0.0 

IF(T.LE. 14500.) GO TO 50 
PB= 0.001 *(T- 14500.) 

PA= 1.0-PB 
50 CONTINUE 
C 

T2= T*T 
T3= T2*T 
T4= T3*T 
TOV= 1.0/T 

C Compute properties for 11 species 
DO 65 1= 1,11 
DO 60 J= 1,6 

60 P(J)= PA*COEF(I,K,J)+ PB*COEF(I,L,J) 

HI(I)= UNIR*T*(P( 1)+ 0.5*P(2)*T+ P(3)<T2/3.+ 0.25*P(4)T3 
1 + 0.2*P(5)n , 4+ P(6)*TOV) 

CPI(I)= UNIR*(P( 1) + P(2)*T+ P(3)*T2+ P(4)T3+ P(5)*T4) 
65 CONTINUE 
RETURN 
END 
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Appendix C 


Electron and Vibrational Heat Fluxes and Associated Thermal Conductivities for 
Flows with Thermal Nonequilibrium 

Electron and vibrational heat fluxes are required in the electron energy and vibra- 
tional energy conservation equations in a multi-temperature gas model (refs. 20 and 
46). The k th component of heat flux vector <7* resulting from the collisions between 
electrons only is given by 

37 ' 

<7*= (CD 


where p, ( =m t n ,) is the density and h, (=(5/2 )kT,/m,) is the enthalpy of electrons, m t 
is the electron mass,! n, is the electron number density, k is Boltzmann’s constant, V* 
is the jfcth component of electron diffusion velocity and K' is the thermal conductivity 
of electrons resulting from collisions between electrons only. In a flowfield analysis, 
the diffusion velocity V k is generally obtained from the electron concentration gradient 
by using the relation 


V k = -D 

v i t] 


3ln(C,) 


(C2) 


where D tj is the binary diffusion coefficient of electrons diffusing into the species j and 
C. is the electron mass fraction. 

The thermal conductivity K, (appearing in equation (Cl)) due to collisions 
between electrons only may be obtained from K, (defined through equation (54)) by 
employing the relation 


K.= K, 


j= 1 


(C3) 


The £th component of vibrational heat flux vector <7* (obtained from molecule- 
molecule collisions only) may be written as 

3 T 

+ E Pi h, vtb V * ( C4) 

uX i = moUcuUs 

where A,.,,* is the vibrational component of enthalpy of the molecules, V* is the £th 
component of diffusion velocity of the molecules, and K’ lb is the thermal conductivity 
due to the collisions between molecules only. 

Referring to Fig. 2, h iiVib for the partially excited vibrational energy mode may be 
obtained from 

= [h,] Tvih - l(C Pti ), + (C^UKT^- 7V,,)- (&h/)T n , (C5) 

In equation (C5), the vibrational component of enthalpy for species i is obtained by 
using the curve-fit relation of equation (14) for the specific enthalpy A, evaluated at 
temperature T nb and subtracting out the contributions from the translational and rota- 
tional enthalpies evaluated at T vlb as well as the enthalpy of formation. Equation (C5) 


t Electron mass m, is related to the electron molecular weight M t through the relation 
m, = M ,/N A , where N A = 6.02x 10 23 molecules / gm — mole is the Avogadro’s number. 
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can be employed up to the characteristic vibrational temperature ©,>* for the species 
i, i.e., for r vil> £ ©i.vib- For higher temperatures (T vi> > ©,„*), the vibrational energy 
mode is fully excited and h iiVib can be obtained from 

= R uni*T vU,! M i (C6) 


The thermal conductivity K'y, (also known as the V -V thermal conductivity) in 
equation (C4) is obtained from the overall vibrational thermal conductivity K vib 
(defined in equation (52)) by considering the contribution of molecule-molecule colli- 
sions only. 


K* = K 


vib 


i s molecules 


S Xj&W(T vib) 

}- molecules 

10 

ZxjA?\T vib) + x t AP(T t) 
/• i 


(C7) 
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Table I. Third Body Efficiences Relative to Argon 


Efficiences relative to argon of 


Catalytic 

Bodies 

Z (y-.V5) # i 

0 2 
i= 1 

yv 2 

2 

N 

3 

0 

4 

NO 

5 

NO* 

6 

01 

7 

Nt 

8 

0 * 

9 

N* 

10 

M, 

l.i 

9 

2 

25 

1 

1 

0 

0 

0 

0 

0 

m 2 

2,i 

1 

2.5 

1 

0 

1 

0 

0 

0 

0 

0 

m 3 

3,i 

1 

1 

20 

20 

20 

0 

0 

0 

0 

0 

m 4 

4,i 

4 

1 

0 

0 

0 

0 

0 

0 

0 

0 

e~ 

5,i 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 
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Table II. Chemical Reactions and Rate Coefficients 


No.f 

r 

Reaction 

Forward rate coefficient, 

kf ,r* 

cm 3 / mole- sec 

Backward rate coefficient, Third body 

*»,r. M 

cm} ! mole- sec or cm 6 ! mole 1 - sec 

1 

0 2 + M x £ 20 + M, 

3 .6 1 x 1 0 18 7 - *°exp (- 5 .94x 1 0 4 / 7 ) 

3.01xl0 15 7-° 5 

0,N.0 2 ,N 2 .N0 

2 

Ni+M 2 £2N + Mi 

1.92xl0 n 7- 05 exp(-1.131xl0 5 /7) 

1.09x l0 16 7-°- 5 

0,02,Ni,N0 

3 

N 2 +N £ 2 N+N 

4.15x 10“r-‘- 5 exp(- 1.3 lx 10 5 /T) 

2.32xl0 21 T- 1 - 5 


4 

NO + MiZtN + O + M) 

3.97x 10 2 °7-' ,J exp(-7.56x 10 4 /7) 

l.OlxlO^T* 15 

o,n,o 2 ,n 2 >no 

5 

NO + O £ 0 2 +jV 

3.18x 10*7 10 exp(- 1.97x 10 4 /7) 

9.63xl0**7 a5 exp(-3.6xl0 3 /7) 


6 

N 2 +0 £ NO+N 

6.75x 10‘ 3 exp(-3.75x 10 4 /7) 

1.5xl0 13 


7 

N + O ;±NO++e- 

9.03x 10 9 7°- 5 exp(-3.24x 10 4 /7) 

1.80x10**7- 10 


8 

O + e' £ 0*+e~+e~ 

(3.6± 1 ,2)x 10 31 7" 191 e*p(- 1.58x 10 5 /7 ) 

(2.2± 0.7)x lO^- 4 - 5 


9 

N + e~ A/*+«-+g- 

(1.1± 0.4)xl0 32 7- JU «p(-1.69xl0 5 /7) 

(2.2± 0.7)xl0 40 7- 4 - 5 


10 

O + O £ 0£ + e" 

(1.6± 0.4)x 10 l7 7'*° ,M ex/>(-8.08x 10 4 /7) 

(8.02± 2.0)xl0 2l 7-*- 5 


11 

0 + 0 2 ^ 0 2 +0* 

2.92x 10**7 - ‘“exp (- 2.8x 10 4 /7 ) 

7.8x1 0**7 05 


12 

(v 2 +(v* ^yv+yvj 

2.02xl0*>7 a8 ‘exp(-1.3xl0 4 /7) 

7.8xl0"7 a5 


13 

/V + iV /V 2 +e~ 

(1.4+ 0.3)xl0"exp(.-(>.lSxlO i /T) 

(1.5± 0.5)x 10 22 7 _ 1-3 


14 

O 2 + /V j £ NO + NO r + e~ 

1.38x io 2 ^-*' 84 exp (1.4 lx 10 5 /7) 

l.Ox 10 24 7" 25 


15 

NO + Mt^NO' + e' + Mt 

2.2x 10 l5 r- OM exp(- 1.08x 1 0 3 / 7* ) 

2.2x 1 0*7 " 2-5 

02 .N 2 

16 

0 + NO* £ NO + O* 

3.63x 10 15 r-°- 6 exp(-5.08xl0 4 /D 

1.5xl0 13 


17 

N 2 + O'' 0 + jV 2 * 

3.4x 10 19 r- 2 -°exp(-2.3x 10 4 /7) 

2.48x 10 i9 7" 22 


18 

N +NO~ ^*NO+N* 

1 .Ox 10 l9 T-°- n exp(- 6.1 x 10 4 /T ) 

4.8x 10 14 


19 

0 2 +NO~ ZNO + O 2 

1.8xl0‘ 5 7°* 7 exp(-3.3xl0 4 /7) 

1.8xl0 13 7 a5 


20 

0 +NO" £ O 2 +N* 

1.34x 10‘ 3 7 a3 ‘exp(-7.727x 10 4 /7) 

l.Ox 10 14 



+ Reactions and reaction rates for Nos. 1-7 are from Blottner’s (ref. 13) 7-species chemical model and Nos. 8-20 are from 
reference 17. 
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Table III. Constants for Polynomial Curve-fits of Thermodynamic Properties 

(300 Ts 30000 K)t 



0.36748£+ 01 -0.12081£-02 0.23240£-05 -0.63218£-09 -0.22577£-12 -0.10612£ + 04 0.23580£+01 

0.28963£ + 01 0.15155£-02 -0.57235£-06 0.99807£-10 -0.65224£-14 -0.90586£+03 0.61615£ + 01 

N 2 0.37270£+ 01 0.46840£-03 -0.11400£-06 0.11540£-10 -0.32930£-15 -0.10430£ + 04 0.12940£+01 

0.96377£+ 01 -0.25728£-02 0.33020£-06 -0.14315£-10 0.20333£-15 -0.10430£+04 -0.37587£+02 

— 0.51 68 1 £-h 01 0.23337£-02 -0.12953£-06 0.27872£-ll -0.21360£-16 -0.10430£ + 04 0.66217£+02 

0.36256£+01 -0.18782£-02 0.705S5E-05 -0.67635E-08 0.21556£-11 -0.10475£ + 04 0.43053£ + 01 

0.36219£+01 0.73618£-03 -0.19652£-06 0.36202£-10 -0.28946£-14 - 0.12020£ + 04 0.36151E+01 

0 2 0.37210£+ 01 0.42540£-03 -0.28350£-07 0.60500£-12 -0.51860£-17 - 0.10440£ + 04 0.32540£+01 

0.34867£+01 0.52384£- 03 -0.39123£-07 0.10094£-11 -0.88718£-17 -0.10440£ + 04 0.48179£ + 01 

0.39620£+ 01 0.39446£- 03 -0.29506£-07 0.73975£-12 -0.64209£-17 -0.10440£ + 04 0.13985£ + 01 

0.25031£+01 -0.21800£-04 0.54205£-07 -0.56476£-10 0.20999£-13 0.56099£ + 05 0.41676£+01 

0.24503£ + 01 0.1 0661 £-03 -0.74653£-07 0.18797£-10 -0.10260£-14 0.56116£+05 0.44488£ + 01 

N 0.27480£ + 01 -0.39090£-03 0.13380£-06 -0.1 1910£— 10 0.33690£-15 0.56090£ + 05 0.28720£+01 

-0.12280£+0l 0.19268£-02 -0.24370£-06 0.12193£-10 -0.19918£-15 0.56090£ + 05 0.28469£ + 02 

0.15520£+02 -0.38858£- 02 0.32288£-06 -0.96053£-U 0.95472£-16 0.56090£+05 -0.88120£ + 02 

0.29464£+01 -0.16382£-02 0.24210£-05 -0.16028£-08 0.38907£-l2 0.29148£ + 05 0.29640£+01 

0.25421£+ 01 - 0.2755 1£- 04 -0.31028£-08 0.45511£-11 -0.43680£-15 0.29231£ + 05 0.49203£ + 01 

O 0.25460£ + 01 -0.59520£-04 0.27010£-07 -0.27980£-11 0.93800£-16 0.29150£ + 05 0.50490£+01 

-0.97871£-02 0.12450£-02 -0.16154£-06 0.80380£-ll -0.12624£-15 0.29150£ + 05 0.21711£+02 

0.16428£+02 -0.39313£-02 0.29840£-06 -0.81613£-11 0.75004£-16 0.29150£+05 -0.94358£+02 

0.40459£ + 01 -0.34182£- 02 0.79819£-05 - 0.61 139£-08 0.15919£- 1 1 0.97454£+04 0.29975£+01 

0.31890£+01 0.13382£- 02 -0.52899£-06 0.95919£— 10 -0.64848£-14 0.98283£ + 04 0.67458£+01 

NO 0.38450£+01 0.25210£-03 -0.26580£-07 0.21620£-11 — 0.63810£- 16 0.97640£ + 04 0.32120£ + 01 

0.43309£+01 -0.58086£-04 0.28059£-07 -0.15694£-11 0.24104£-16 0.97640£ + 04 0.10735£ + 00 

0.23507£ + 01 0.58643£-03 -0.31316£-07 0.60495£-12 -0.40557£-17 0.97640£+04 0.14026£ + 02 

0.36685£+01 -0.11545£-02 0.21756£-05 -0.48227£-09 -0.27848£-12 0.11803£ + 06 0.31779£+01 

0.28886£+ 01 0.15217£-02 -0.57531£-06 0.10051£-09 -0.66044£-14 0.11819£ + 06 0.70027£ + 01 

NO* 0.22142£+ 01 0.17761ZT- 02 -0.43039£-06 0.41738£- 10 -0.12829£-14 0.11819£ + 06 0.11268£ + 02 

-0.33240£ + 01 0.24420E- 02 -0.19057£-06 0.68580£- 1 1 -0.99112£-16 0.11819£ + 06 0.51864£+02 

-0.43488£ + 01 0.24012£-02 -0.14460£-06 0.33813£-11 -0.28255£-16 0.11819£ + 06 0.60082£ + 02 

0.25000£ + 01 0.00000£ + 00 0.00000£ + 00 0.00000£ + 00 0.00000£ + 00 -0.74538£ + 03 -0.11734£ + 02 

0.25000£+01 0.00000£ + 00 0.00000£ + 00 0.00000£+00 0.00000£ + 00 -0.74538£ + 03 -0.11734£ + 02 

e* 0.25080£+ 01 -0.63320£-05 0.13640£-08 -0.10940£-12 0.29340£-17 - 0.74500£ + 03 -0.12080£ + 02 

0.25001£+ 01 -0.31128E-09 0.35721£- 13 -0.16037£-17 0.25071£-22 - 0.74500£ + 03 -0.11734£ + 02 

0.25001£+ 01 0.30158£- 09 -0.22620£- 13 0.66734£-18 -0.68917£-23 -0.74500£ + 03 -0.11734£ + 02 

0.27270£+01 -0.28200£-03 0.1 1050£- 06 -0.15510£- 10 0.78470£-15 0.22540£ + 06 0.36450£+01 

0.27270£+01 -0.28200E-03 0.11050£-06 -0.15510£- 10 0.78470£-15 0.22540£ + 06 0.36450£ + 01 

N* 0.24990£+01 -0.37250£-05 0.11470£-07 -0.1 1020£- 1 1 0.30780£-16 0.22540£ + 06 0.49S00£+01 

0.23856£+ 01 0.83495£-04 -0.58815£-08 0.18850£- 12 -0.16120£-17 0.22540£ + 06 0.56462£+01 

0.22286£ + 01 0.12458£-03 -0.87636£-08 0.26204£-12 -0.21674£— 17 0.22540£ + 06 0.67811£ + 01 
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0 + 


0.24985E + 01 
0.25061£+01 
0.29440E+ 01 
0.12784E+ 01 
0.12889E+ 01 


0.1 141 l£-04 
-0.14464E-04 
-0.41080E-03 
0.40866£-03 
0.43343£-03 


-0.2976l£-07 

0.12446£-07 

0.91560£-07 

-0.21731£-07 

-0.26758£-07 


Q.32247E- 10 
-0.46858£ - 11 
-0.58480E- 1 1 
0.33252E- 12 
0.62159E- 12 


-0.12376£ - 13 
0.65549£- 15 
0.1 1900£- 15 
0.63 160£- 18 
-0.45131£- 17 


0.18795£ + 06 
0.18795£ + 06 
0.18790E + 06 
0.18790E+06 
0.18790£ + 06 


0.43864£ + 01 
0.43480E + 01 
0.17500E + 01 
0.12761E + 02 
0.12604E + 02 


Ni 


0.33970E+ 01 
0.33970E + 01 
0.33780£+ 01 
0.43942E + 01 
0.39493£ + 01 


0.45250£- 03 
0.45250£- 03 
0.86290£-03 
0.18868E- 03 
0.36795£-03 


0.12720E-06 

0.12720E-06 

-0.12760£-06 

- 0.7 1272£ - 08 

- 0.269 10£- 07 


-0.38790E-10 
-0.38790E- 10 
0.80870E - 11 
-0.17511E-12 
0.67110E- 12 


0.24590£- 14 
0.24590E- 14 
-0.18800£ - 15 
0.67176£- 17 
-0.58244£-17 


0.18260£ + 06 
0.18260£+06 
0.18260£ + 06 
0.18260£ + 06 
0.18260£+06 


0.42050£ + 01 
0.42050£ + 01 
0.40730E + 01 
-0.23693E+01 
0.65472E + 00 


o i 


0.32430£+01 
Q.32430E + 01 
0.51690£+ 01 
-0.28G17E+ 00 
0.20445E + 01 


0.1 1740£- 02 
0.1 1740E-02 
-0.86200£-03 
0.16674E-02 
0.10313£- 02 


-0.39000E-06 
-0.39000E-06 
0.204 10£- 06 
-0.12107£-06 
-0.74046£-07 


0.54370E-10 
0.54370E - 10 
-0.13000E- 10 
0.321 13£- 11 
0.19257E - 11 


-0.23920£- 14 
- 0.23920£ - 1 4 
0.24940E- 15 
-0.28349E - 16 
-0.17461£- 16 


0.14000£ + 06 
0.14000£ + 06 
0.14000£ + 06 
0.14000£ + 06 
0.14000£ + 06 


0.59250E + 01 
0.59250E + 01 
-0.52960£ + 01 
0.31013£ + 02 
0.143 10£ + 02 


t There are five rows of constants provided for each species which correspond to the following five temperature ranges, 
respectively. 


300 K < T £ 1000 K 
1000 K £ T £ 6000 K 
6000 K < T < 15000 K 
15000 K £ T < 25000 K 
25000 K < T <; 30000 K 

For temperatures less than 300 K. (which may arise in the freestream for certain flight trajectories (refs. 15,52)), the 
specific heats of O >Oi % N , and TV 2 are practically constant. For such cases, the specific enthalpy may be obtained from the 
relation: 


h;= C,j(T -T n ,)+^hf) Tnf 

with T~f = 298 K and C p 0 = 5.44 call gm- mole- K , C p0 2 = 7.00 call gm- mole- K , C pN = 4.97 call gm-mole- K , 
C P 'H 1 = 6.96 cal! gm- mole- K . 
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Table IV. Constants for Viscosityf Curve-Fits 
(1000 K < T £ 30000 K)tt 


Speciesl -4^. B |i . D ^ E^. 


N 2 

0.0 

02 

0.0 

N 

0.0 

O 

0.0 

NO 

0.0 

NO + 

0.0913 

e~ 

0.0899 

N + 

0.0895 

0 + 

0.0912 


0.0897 

0} 

0.0908 


0.0 

0.0203 

0.0 

0.0484 

0.0 

0.0120 

0.0 

0.0205 

0.0 

0.0452 

-3.3178 

45.1426 

-3.2731 

44.5782 

-3.2573 

44.3511 

-3.3154 

45.1290 

-3.2618 

44.4079 

-3.3020 

44.9511 


0.4329 

-11.8153 

-0.1455 

-8.9231 

0.5930 

-12.3805 

0.4257 

-11.5803 

-0.0609 

-9.4596 

-270.3948 

586.3300 

-267.2522 

574.4149 

-265.8276 

576.1313 

-270.4211 

586.2903 

-266.1462 

577.1449 

-269.3877 

584.4130 


t Viscosity is obtained in gm! cm- sec 

tt For temperatures less than 1000 K, Sutherland’s viscosity law for air may be used for all 
species. 

$ The ionic species viscosities given here are for an electron pressure of 1 atm. For different 
electron pressures, these values should be corrected by the formula given in equation (24a) in 
the main text. 
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Table V. Constants for Frozen Thermal Conductivity! Curve-Fits 
(1000 K< T< 30000 K)!* 


Species* A K/ . B K[ . C Kf . D Kf . E Kf . 


n 2 

0.0418 

0 2 

0.0776 

N 

0.0 

0 

0.0 

NO 

0.0263 

NO + 

0.0 

e~ 

0.0908 

N + 

0.1023 

0 + 

0.0534 

*2 

0.0918 

O z 

-0.0125 


-1.2720 

14.4571 

-2.5007 

30.1390 

0.0 

0.0169 

0.0 

0.0319 

-0.8130 

9.4203 

0.0431 

-0.5477 

-3.3046 

44.9877 

-3.6925 

49.8997 

-2.0710 

29.9141 

-3.3760 

46.2945 

0.1276 

2.5428 


-71.9660 

122.5745 

-160.1758 

307.3634 

0.5373 

-12.8682 

0.2485 

-11.6657 

-47.6841 

79.4139 

2.1912 

-17.8610 

-269.6002 

588.9412 

-297.1291 

641.3168 

-188.4951 

421.3914 

-278.4425 

604.2017 

-37.7507 

111.8110 


t Thermal conductivity is obtained in call cm- sec-K 

ft For temperatures lower than 1000 K, Sutherland’s law for thermal conductivity of air may 
be used for each species. 

* The ionic species frozen thermal conductivities given here are for an electron pressure of 1 
atm. For different electron pressure, these values should be corrected by employing the for- 
mula given in equation (24a) in the main text. 
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Table VI. Constants for Diffusion Coefficientf Curve-Fits 


No. ft 

Interaction 
Pair (i-j) 

A °>i 

B *i 


% 

Temperature 
Range, K| 

1 

n 2 - n 2 

0.0 

0.0112 

1.6182 

-11.3091 


2 

0 2 -N 2 

0.0 

0.0465 

0.9271 

-8.1137 


3 

O 2- 

0.0 

0.0410 

1.0023 

-8.3597 


4 

N - N 2 

0.0 

0.0195 

1.4880 

-10.3654 


5 

N - 0 2 

0.0 

0.0179 

1.4848 

-10.2810 


6 

N - N 

0.0 

0.0033 

1.5572 

-11.1616 


7 

o-n 2 

0.0 

0.0140 

1.5824 

-10.8819 


8 

0-0 2 

0.0 

0.0226 

1.3700 

-9.6631 


9 

0 - N 

0.0 

-0.0048 

1.9195 

-11.9261 


10 

0-0 

0.0 

0.0034 

1.5572 

-11.1729 


11 

no-n 2 

0.0 

0.0291 

1.2676 

-9.6878 


12 

NO - 0 2 

0.0 

0.0438 

0.9647 

-8.2380 


13 

NO -N 

0.0 

0.0185 

1.4882 

-10.3301 


14 

NO -O 

0.0 

0.0179 

1.4848 

-10.3155 


15 

NO-NO 

0.0 

0.0364 

1.1176 

-8.9695 


16 

no*-n 2 

0.0 

0.0 

1.9000 

-13.3343 


17 

NO* - 02 

0.0 

0.0 

1.9001 

-13.3677 


18 

NO* - N 

0.0 

0.0 

1.8999 

-13.1254 


19 

NO* -O 

0.0 

0.0 

1.9000 

-13.1701 


20 

NO* -NO 

0.0 

0.0047 

1.5552 

-11.3713 


21 

NO* -NO' 

-0.1251 

3.5135 

-29.7280 

74.1550 


22 

e~-N 2 

-0.1147 

2.8945 

-23.0085 

65.9815 


23 

e~ - 0 2 

-0.0241 

0.3464 

0.1136 

-1.3848 

1000-9000 



-0.0029 

0.0856 

0.6655 

-0.8205 

9000-30000 

24 

e~ - N 

0.0 

0.0 

1.5000 

-2.9987 


25 

e~ - O 

0.0581 

-1.5975 

15.4508 

-40.7370 


26 

e ~-NO 

0.2202 

-5.2261 

42.0630 

-106.0937 

1000-8000 



0.2871 

-8.3759 

82.8802 

-267.0227 

8000-30000 

27 

e '-NO* 

-0.1251 

3.5134 

-29.7272 

79.2610 


28 

e~ - e~ 

-0.1251 

3.5136 

-29.7290 

79.6126 


29 

N*- N 2 

0.0 

0.0 

1.9000 

-13.1144 


30 

N*- O j 

0.0 

0.0 

1.9000 

-13.1357 


31 

N*- N 

0.0 

0.0033 

1.5572 

-11.1616 


32 

N*- 0 

0.0 

0.0 

1.9000 

-13.0028 


33 

N* - NO 

0.0 

0.0 

1.8999 

-13.1254 


34 

N * - -NO* 

-0.1251 

3.5135 

-29.7285 

74.3825 


35 

N* - e~ 

-0.1251 

3.5134 

-29.7272 

79.2611 


36 

N*-N* 

-0.1251 

3.5134 

-29.7274 

74.5342 


37 

o*-n 2 

0.0 

0.0 

1.9000 

-13.1578 


38 

O* -0 2 

0.0 

0.0 

1.9000 

-13.1810 


39 

0 * -N 

0.0 

0.0 

1.9000 

-13.0028 


40 

O* -O 

0.0 

0.0034 

1.5572 

-11.1729 


41 

O*- NO 

0.0 

0.0 

1.9000 

-13.1701 


42 

O*- NO* 

-0.1251 

3.5133 

-29.7268 

74.3330 


43 

0*-e~ 

-0.1251 

3.5134 

-29.7274 

79.2616 
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44 

0 + 

- N + 

- 0.1251 

3.5134 

- 29.7275 

74.5023 

45 

0 + 

- 0 + 

- 0.1251 

3.5134 

- 29.7277 

74.4687 

46 

Nt 

-n 2 

0.0 

0.0 

1.9000 

- 13.3173 

47 

Nl 

- O 2 

0.0 

0.0 

1.9000 

- 13.3495 

48 

N 2 

- N 

0.0 

0.0 

1.9000 

- 13.1144 

49 

n; 

- O 

0.0 

0.0 

1.9000 

- 13.1578 

50 

N 2 + 

-NO 

0.0 

0.0 

1.9000 

- 13.3343 

51 

*2 

- NO + 

- 0.1251 

3.5134 

- 29.7279 

74.1721 

52 

N 2 

- e~ 

- 0.1251 

3.5134 

- 29.7273 

79.2613 

53 

N 2 

- N + 

- 0.1251 

3.5135 

- 29.7288 

74.3947 

54 

n; 

- 0 + 

- 0.1251 

3.5133 

- 29.7269 

74.3453 

55 

N 2 

-Ni 

- 0.1251 

3.5135 

- 29.7282 

74.1899 

56 

oi 

— N 2 

0.0 

0.0 

1.9000 

- 13.3173 

57 

oi 

— O 2 

0.0 

0.0 

1.9000 

- 13.3495 

58 

oi 

-N~ 

0.0 

0.0 

1.9000 

- 13.1144 

59 

O 2 + 

- O 

0.0 

0.0 

1.9000 

- 13.1578 

60 

oi 

-NO 

0.0 

0.0 

1.9000 

- 13.3343 

61 

oi 

- NO + 

- 0.1251 

3.5134 

- 29.7279 

74.1721 

62 

oi 

- e~ 

- 0.1251 

3.5134 

- 29.7273 

79.2613 

63 

oi 

- N + 

- 0.1251 

3.5135 

- 29.7288 

74.3947 

64 

oi 

- 0 + 

- 0.1251 

3.5133 

- 29.7269 

74.3453 

65 

oi 

-Ni 

- 0.1251 

3.5135 

- 29.7282 

74.1899 

66 

oi 

-oi 

- 0.1251 

3.5135 

- 29.7282 

74.1899 


t Diffusion coefficients are obtained in cm 2 - atm! sec. Diffusion coefficients obtained from these curve- 
fits are for an electron pressure of 1 atm. For different electron pressures, the cross sections should be 
corrected by the formula given in the main text when the interacting pair of species are both ions or elec- 
trons or a combination of the two. Note that diffusion coefficients for N { and Oi are taken to be the 
same. 

ft Cross section Nos. 1-15 are used in a 5-species air model and Nos. 1-28 in a 7-species model, 
t The temperature range for all curve fits is 1000 < T £ 30000 K, except where noted. 
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Table VII. Curve-Fit Constants for Collision Cross-Section, 


No.tt 

Interaction 
Pair (i-j) 



w V 


Temperature 
Range, K$ 

1 

n 2 -n 2 

0.0 

-0.0112 

-0.1182 

4.8464 


2 

0 2 “ N 2 

0.0 

-0.0465 

0.5729 

1.6185 


3 

0 2 — 0 2 

0.0 

-0.0410 

0.4977 

1.8302 


4 

N - N 2 

0.0 

-0.0194 

0.0119 

4.1055 


5 

J V - 0 2 

0.0 

-0.0179 

0.0152 

3.9996 


6 

N - N 

0.0 

-0.0033 

-0.0572 

5.0452 


7 

o-n 2 

0.0 

-0.0139 

-0.0825 

4.5785 


8 

0-0 2 

0.0 

-0.0226 

0.1300 

3.3363 


9 

O-N 

0.0 

0.0048 

-0.4195 

5.7774 


10 

0-0 

0.0 

-0.0034 

-0.0572 

4.9901 


11 

no-n 2 

0.0 

-0.0291 

0.2324 

3.2082 


12 

NO - 0 2 

0.0 

-0.0438 

0.5352 

1.7252 


13 

NO -N 

0.0 

-0.0185 

0.0118 

4.0590 


14 

NO -0 

0.0 

-0.0179 

0.0152 

3.9996 


15 

NO - NO 

0.0 

-0.0364 

0.3825 

2.4718 


16 

no*-n 2 

0.0 

0.0 

- 0.4000 

6.8543 


17 

no*-o 2 

0.0 

0.0 

-0.4000 

6.8543 


18 

NO* ~n' 

0.0 

0.0 

-0.4000 

6.8543 


19 

NO* -0 

0.0 

0.0 

-0.4000 

6.8543 


20 

NO* -NO 

0.0 

-0.0047 

-0.0551 

4.8737 


21 

NO* -NO* 

0.1251 

-3.5134 

31.2277 

-80.6515 


22 

e~ -N 2 

0.1147 

-2.8945 

24.5080 

-67.3691 


23 

- 0 2 

0.0241 

-0.3467 

1.3887 

-0.0110 

1000-9000 



0.0025 

-0.0742 

0.7235 

-0.2116 

9000-30000 

24 

e~-N 

0.0 

0.0 

0.0 

1.6094 


25 

e~ - 0 

0.0164 

-0.2431 

1.1231 

-1.5561 

1000-9000 



-0.2027 

5.6428 

-51.5646 

155.6091 

9000-30000 

26 

e~ - NO 

- 0.2202 

5.2265 

-40.5659 

104.7126 

1000-8000 



-0.2871 

8.3757 

-81.3787 

265.6292 

8000-30000 

27 

e ~-NO* 

0.1251 

-3.5134 

31.2277 

-80.6515 


28 

e~ - e~ 

0.1251 

-3.5134 

31.2277 

-80.6515 


29 

N*- N 2 

0.0 

0.0 

- 0.4000 

6.8543 


30 

n*-o 2 

0.0 

0.0 

-0.4000 

6.8543 


31 

N* - N 

0.0 

-0.0033 

-0.0572 

5.0452 


32 

N*-0 

0.0 

0.0 

- 0.4000 

6.8543 


33 

N*- -NO 

0.0 

0.0 

- 0.4000 

6.8543 


34 

N * -NO* 

0.1251 

-3.5134 

31.2277 

- 80.6515 


35 

N*- e~ 

0.1251 

-3.5134 

31.2277 

-80.6515 


36 

N*-N* 

0.1251 

-3.5134 

31.2277 

-80.6515 


37 

o*-n 2 

0.0 

0.0 

-0.4000 

6.8543 


38 

0*- 0 2 

0.0 

0.0 

-0.4000 

6.8543 


39 

0* - N 

0.0 

0.0 

-0.4000 

6.8543 


40 

0*-0 

0.0 

-0.0034 

-0.0572 

4.9901 


41 

0* - NO 

0.0 

0.0 

-0.4000 

6.8543 


42 

0* - NO* 

0.1251 

-3.5134 

31.2277 

-80.6515 
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43 

O* 

- e 

0.1251 

- 3.5134 

31.2277 

- 80.6515 

44 

O* 

-N + 

0.1251 

- 3.5134 

31.2277 

- 80.6515 

45 

0 + 

- 0 + 

0.1251 

- 3.5134 

31.2277 

- 80.6515 

46 

N 2 

-n 2 

0.0 

0.0 

- 0.4000 

6.8543 

47 

N 2 + 

- 0 2 

0.0 

0.0 

- 0.4000 

6.8543 

48 

*2 

- N 

0.0 

0.0 

- 0.4000 

6.8543 

49 

N 2 + 

- O 

0.0 

0.0 

- 0.4000 

6.8543 

50 

N 2 

- NO 

0.0 

0.0 

- 0.4000 

6.8543 

51 

n; 

- NO + 

0.1251 

- 3.5134 

31.2277 

- 80.6515 

52 

N 2 + 

- e~ 

0.1251 

- 3.5134 

31.2277 

- 80.6515 

53 

Nt 

- N + 

0.1251 

- 3.5134 

31.2277 

- 80.6515 

54 

*2 

- 0 + 

0.1251 

- 3.5134 

31.2277 

- 80.6515 

55 


2 

0.1251 

- 3.5134 

31.2277 

- 80.6515 

56 

Of 

~ N 2 

0.0 

0.0 

- 0.4000 

6.8543 

57 

o 2 + 

- o 2 

0.0 

0.0 

- 0.4000 

6.8543 

58 

o 2 + 

-N 

0.0 

0.0 

- 0.4000 

6.8543 

59 

ot 

- O 

0.0 

0.0 

- 0.4000 

6.8543 

60 

Of 

-NO 

0.0 

0.0 

- 0.4000 

6.8543 

61 

°2 

-NO* 

0.1251 

- 3.5134 

31.2277 

- 80.6515 

62 

°2 

- e~ 

0.1251 

- 3.5134 

31.2277 

- 80.6515 

63 

0$ 

- N* 

0.1251 

- 3.5134 

31.2277 

- 80.6515 

64 

ot 

- O* 

0.1251 

- 3.5134 

31.2277 

- 80.6515 

65 

0 2 + 

-N± 

0.1251 

- 3.5134 

31.2277 

- 80.6515 

66 

o; 

- O! 

0.1251 

- 3.5134 

31.2277 

- 80.6515 


t Cross sections are obtained in A 2 ; 1 A 2 = 10" 16 cm 2 . Collision cross sections obtained from these 
curve-fits are for an electron pressure of 1 aim. For different electron pressures, the cross sections should 
be corrected by the formula given in the main text when the interacting pair of species are both ions or 
electrons or a combination of the two. Note that cross sections for N { and 01 are taken to be the same. 

ft Cross section Nos. 1-15 are used in a 5-species air model and Nos. 1-28 in a 7-species model. 

t The temperature range for all curve fits is 1000 < T £ 30000 K, except where noted. 
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Table VIII. Curve-Fit Constants for Collision Cross-Section, f 


No. ft 

Interaction 
Pair (i-j) 

A af 2) 

B n^ 2 > 

c sr?. 2 ' 2 ) 

O 77(2,2) 
4 ■; 

Temperature 
Range, Kf 

1 

n 2 - n 2 

0.0 

-0.0203 

0.0683 

4.0900 


2 

o 2 - n 2 

0.0 

-0.0558 

0.7590 

0.8955 


3 

0 2 - 0 2 

0.0 

-0.0485 

0.6475 

1.2607 


4 

N - N 2 

0.0 

-0.0190 

0.0239 

4.1782 


5 

N - 0 2 

0.0 

-0.0203 

0.0730 

3.8818 


6 

N -N 

0.0 

-0.0118 

-0.0960 

4.3252 


7 

o - n 2 

0.0 

-0.0169 

-0.0143 

4.4195 


8 

0-0 2 

0.0 

-0.0247 

0.1783 

3.2517 


9 

0 - N 

0.0 

0.0065 

-0.4467 

6.0426 


10 

0-0 

0.0 

-0.0207 

0.0780 

3.5658 


11 

no-n 2 

0.0 

-0.0385 

0.4226 

2.4507 


12 

NO - 0 2 

0.0 

-0.0522 

0.7045 

1.0738 


13 

NO -N 

0.0 

-0.0196 

0.0478 

4.0321 


14 

NO - 0 

0.0 

-0.0203 

0.0730 

3.8818 


15 

NO ~ NO 

0.0 

-0.0453 

0.5624 

1.7669 


16 

no*-n 2 

0.0 

0.0 

-0.4000 

6.7760 


17 

no* -o 2 

0.0 

0.0 

-0.4000 

6.7760 


18 

NO * - N 

0.0 

0.0 

-0.4000 

6.7760 


19 

NO * - 0 

0.0 

0.0 

-0.4000 

6.7760 


20 

NO* -NO 

0.0 

0.0 

-0.4000 

6.7760 


21 

NO* -NO* 

0.1251 

-3.5135 

31.2281 

-80.1163 


22 

e~-N 2 

0.1147 

-2.8945 

24.5080 

-67.3691 


23 

e ~ - 0 2 

0.0241 

-0.3467 

1.3887 

-0.0110 

1000-9000 



0.0025 

-0.0742 

0.7235 

-0.2116 

9000-30000 

24 

e~-N 

0.0 

0.0 

0.0 

1.6094 


25 

e'-0 

0.0164 

-0.2431 

1.1231 

-1.5561 

1000-9000 



-0.2027 

5.6428 

-51.5646 

155.6091 

9000-30000 

26 

e~- NO 

-0.2202 

5.2265 

-40.5659 

104.7126 

1000-8000 



-0.2871 

8.3757 

-81.3787 

265.6292 

8000-30000 

27 

e~ - NO* 

0.1251 

-3.5134 

31.2274 

-80.1684 


28 

e ~ - e ~ 

0.1251 

-3.5134 

31.2274 

-80.1684 


29 

n*-n 2 

0.0 

0.0 

-0.4000 

6.7760 


30 

n*-o 2 

0.0 

0.0 

-0.4000 

6.7760 


31 

N*- N 

0.0 

0.0 

-0.4146 

6.9078 


32 

N*- 0 

0.0 

0.0 

-0.4000 

6.7760 


33 

N*- -NO 

0.0 

0.0 

-0.4000 

6.7760 


34 

N*-NO* 

0.1251 

-3.5135 

31.2281 

-80.1163 


35 

N *- e ~ 

0.1251 

-3.5134 

31.2274 

-80.1684 


36 

N*- N* 

0.1251 

-3.5135 

31.2281 

-80.1163 


37 

o*-n 2 

0.0 

0.0 

-0.4000 

6.7760 


38 

o*-o 2 

0.0 

0.0 

-0.4000 

6.7760 


39 

0* - N 

0.0 

0.0 

-0.4000 

6.7760 


40 

0 * -0 

0.0 

0.0 

-0.4235 

6.7787 


41 

0* -NO 

0.0 

0.0 

-0.4000 

6.7760 


42 

0* - NO* 

0.1251 

-3.5135 

31.2281 

-80.1163 
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43 

O' 

- e~ 

0.1251 

- 3.5134 

31.2274 

- 80.1684 

44 

0 + 

- N + 

0.1251 

- 3.5135 

31.2281 

- 80.1163 

45 

0 + 

- 0 + 

0.1251 

- 3.5135 

31.2281 

- 80.1163 

46 

N Z 

-n 2 

0.0 

0.0 

- 0.4000 

6.7760 

47 

NZ 

-O 2 

0.0 

0.0 

- 0.4000 

6.7760 

48 

N* 

- N 

0.0 

0.0 

- 0.4000 

6.7760 

49 

N 2 

- 0 

0.0 

0.0 

- 0.4000 

6.7760 

50 

*2 

-NO 

0.0 

0.0 

- 0.4000 

6.7760 

51 

*2 

-NO + 

0.1251 

- 3.5135 

31.2281 

- 80.1163 

52 

*2 

- e~ 

0.1251 

- 3.5134 

31.2274 

- 80.1684 

53 

N 2 

- N + 

0.1251 

- 3.5135 

31.2281 

- 80.1163 

54 

nz 

- 0 + 

0.1251 

- 3.5135 

31.2281 

- 80.1163 

55 

NZ 

-NZ 

0.1251 

- 3.5135 

31.2281 

- 80.1163 

56 

oz 

-N 2 

0.0 

0.0 

- 0.4000 

6.7760 

57 

O Z 

- O 2 

0.0 

0.0 

- 0.4000 

6.7760 

58 

oz 

- N 

0.0 

0.0 

- 0.4000 

6.7760 

59 

Ol 

- O 

0.0 

0.0 

- 0.4000 

6.7760 

60 

OZ 

- NO 

0.0 

0.0 

- 0.4000 

6.7760 

61 

oz 

-NO + 

0.1251 

- 3.5135 

31.2281 

- 80.1163 

62 

OZ 

- e~ 

0.1251 

- 3.5135 

31.2281 

- 80.1163 

63 

OZ 

- N + 

0.1251 

- 3.5135 

31.2281 

- 80.1163 

64 

oz 

- 0 ¥ 

0.1251 

- 3.5135 

31.2281 

- 80.1163 

65 

oz 

-NZ 

0.1251 

- 3.5135 

31.2281 

- 80.1163 

66 

OZ 

-oz 

0.1251 

- 3.5135 

31.2281 

- 80.1163 


t Cross sections are obtained in A 2 ; 1 A 2 = 10' 16 cm 2 . Collision cross sections obtained from these 
curve-fits are for an electron pressure of 1 atm. For different electron pressures, the cross sections should 
be corrected by the formula given in the main text when the interacting pair of species are both ions or 
electrons or a combination of the two. Note that cross sections for and O { are taken to be the same. 

ft Cross section xNos. 1-15 are used in a 5-species air model and Nos. 1-28 in a 7-species model. 

X The temperature range for all curve fits is 1000 < T < 30000 K, except where noted. 
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Table IX. Curve-Fit Constants for Collision Cross-Section Ratio, B \ } t 

( 1000 K< T < 30000 K) 


No.tt 

Interaction 
Pair (i-j) 

a r' 

°<i 


c fl ; * 

D b- 

B a 

1 

n 2 - n 2 

0.0 

-0.0073 

0.1444 

-0.5625 

2 

0 2 ~ N 2 

0.0 

-0.0019 

0.0602 

-0.2175 

3 

0 2 ~ 0 2 

0.0 

0.0001 

0.0181 

-0.0306 

4 

N - N 2 

0.0 

0.0043 

-0.0494 

0.2850 

5 

n - o 2 

0.0 

0.0033 

-0.0366 

0.2332 

6 

N -N 

0.0 

0.0002 

0.0002 

0.0537 

7 

o-n 2 

0.0 

0.0042 

-0.0471 

0.2747 

8 

o-o 2 

0.0 

0.0024 

-0.0245 

0.1808 

9 

0 -N 

0.0 

0.0147 

-0.2628 

1.2943 

10 

0-0 

0.0 

0.0002 

0.0 

0.0549 

11 

NO - N 2 

0.0 

-0.0045 

0.1010 

-0.3872 

12 

NO - 0 2 

0.0 

-0.0010 

0.0410 

-0.1312 

13 

NO -N 

0.0 

0.0038 

-0.0425 

0.2574 

14 

NO -O 

0.0 

0.0033 

-0.0366 

0.2332 

15 

NO - NO 

0.0 

-0.0027 

0.0700 

-0.2553 

16 

no * -n 2 

0.0 

0.0 

0.0 

0.1933 

17 

no * - o 2 

0.0 

0.0 

0.0 

0.1933 

18 

NO*-N 

0.0 

0.0 

0.0 

0.1933 

19 

NO* - O 

0.0 

0.0 

0.0 

0.1933 

20 

NO* -NO 

0.0 

0.0003 

-0.0006 

0.0632 

21 

NO* -NO* 

0.0 

0.0 

0.0 

0.4463 

22 

e~ - N 2 

0.0 

0.0 

0.0 

0.0 

23 

e~ - O 2 

0.0 

0.0 

0.0 

0.0 

24 

e~-N 

0.0 

0.0 

0.0 

0.0 

25 

e~ - O 

0.0 

0.0 

0.0 

0.0 

26 

e~ - NO 

0.0 

0.0 

0.0 

0.0 

27 

e~ -NO* 

0.0 

0.0 

0.0 

0.4463 

28 

e~ - e~ 

0.0 

0.0 

0.0 

0.4463 

29 

n *- n 2 

0.0 

0.0 

0.0 

0.1933 

30 

n *- o 2 

0.0 

0.0 

0.0 

0.1933 

31 

N* - N 

0.0 

0.0002 

0.0002 

0.0537 

32 

N* - 0 

0.0 

0.0 

0.0 

0.1933 

33 

N* - NO 

0.0 

0.0 

0.0 

0.1933 

34 

N* -NO* 

0.0 

0.0 

0.0 

0.4463 

35 

N*-e~ 

0.0 

0.0 

0.0 

0.4463 

36 

AT - AT 

0.0 

0.0 

0.0 

0.4463 

37 

0*-N 2 

0.0 

0.0 

0.0 

0.1933 

38 

0~ - 0 2 

0.0 

0.0 

0.0 

0.1933 

39 

0* - N 

0.0 

0.0 

0.0 

0.1933 

40 

0* -0 

0.0 

0.0002 

0.0 

0.0549 

41 

O' -NO 

0.0 

0.0 

0.0 

0.1933 

42 

0* -NO^ 

0.0 

0.0 

0.0 

0.4463 

43 

0*-e' 

0.0 

0.0 

0.0 

0.4463 

44 

0* - N* 

0.0 

0.0 

0.0 

0.4463 
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45 

0 + -0 + 

0.0 

0.0 

0.0 

0.4463 

46 

N 1 - N 2 

0.0 

0.0 

0.0 

0.1933 

47 

ni - 0 2 

0.0 

0.0 

0.0 

0.1933 

48 

Nl - n ' 

0.0 

0.0 

0.0 

0.1933 

49 

N 2 - 0 

0.0 

0.0 

0.0 

0.1933 

50 

N 2 - NO 

0.0 

0.0 

0.0 

0.1933 

51 

N j - NO + 

0.0 

0.0 

0.0 

0.4463 

52 

N 2 - e ~ 

0.0 

0.0 

0.0 

0.4463 

53 

- N + 

0.0 

0.0 

0.0 

0.4463 

54 

Nl - o + 

0.0 


0.0 

0.4463 

55 

N 2 - N $ 

0.0 

0.0 

0.0 

0.4463 

56 

Ot - N 2 

0.0 

0.0 

0.0 

0.1933 

57 

Ot - 0 2 

0.0 

0.0 

0.0 

0.1933 

58 

Ot -N 

0.0 

0.0 

0.0 

0.1933 

59 

oi - 0 

0.0 

0.0 

0.0 

0.1933 

60 

0 2 + - NO 

0.0 

0.0 

0.0 

0.1933 

61 

Of - NO + 

0.0 

0.0 

0.0 

0.4463 

62 

0 2 + - e ~ 

0.0 

0.0 

0.0 

0.4463 

63 

0 2 - N + 

0.0 

0.0 

0.0 

0.4463 

64 

0 2 - 0 + 

0.0 

0.0 

0.0 

0.4463 

65 

Ot -Nl 


0.0 

0.0 

0.4463 

66 

01 - 01 

0.0 

0.0 

0.0 

0.4463 


t The collision cross section ratios are dimensionless parameters and are valid as given for all electron 
pressures. 

ft Cross section Nos. 1-15 are used in a 5-species air model and Nos. 1-28 in a 7-species model. 





Figure 1. Flight stagnation region air chemistry (adapted from ref. 5). 
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Curve fit to McBride-Browne data 
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Figure 3. Curve fit to the specific heat values obtained by McBride et al. (ref. 

26) and Browne (ref. 28) and comparison with Hansen’s values (ref. 
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Temperature, 



Curve fit to McBride-Browne data 



Figure 3. Concluded. 
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Figure 4. Curve fit to the viscosity values obtained by employing the collision 
cross sections based on the data of Yos (ref. 32). 
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(c) Ionized nitric oxide. 
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Figure 5. Curve fit to the frozen thermal conductivity values obtained by 
employing the collision cross sections based on the data of Yos (ref. 



Computed 


-58- 



Figure 5. Continued. 
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Figure 6. Comparison of the viscosity of equilibrium air at 1 atmosphere. 
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Figure 7. Comparison of the frozen thermal conductivity of equilibrium air at 
1 atmosphere. 
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(b) Ion-neutral molecular interaction: NO + <-» £>2 



Computed: 


- 65 - 


O 

O ■ <D 

a. cm « 

o 


CM 

C\T , CM 

o" 0<t 


(c) Ion-ion molecular interaction: <-» NO 
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(d) Electron-ionized-molecule interaction: NO 
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